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ABSTRACT 


Heat  transfer  tests  were  carried  out  using  a  small  enhanced  tube  bundle  in  a  pool  of 
R-113.  By  accurately  instrumenting  five  tubes  within  the  bundle,  both  the  convective  and 
nucleate  boiling  regions  were  studied  in  detail,  with  emphasis  on  the  'bundle'  effect  (i.e. 
the  effect  of  lower  tubes  in  operation  on  upper  tubes  within  the  bundle).  In  addition,  the 
effect  of  surface  history  and  pool  height  on  nucleation  site  activation/deactivation  was 
studied  to  see  how  this  affects  the  overall  heat  transfer  and  in  particular,  the  shape  of  the 
hysteresis  loop.  From  the  results,  recommendations  can  be  made  to  improve  start-up 
procedures  on  shipboard  AC  systems. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  September  1987,  due  to  the  pressing  environmental  issues  of  ozone  depletion  and 
global  warming,  an  international  conference  was  held  in  Montreal,  Canada.  This 
conference  was  viewed  by  many  as  a  major  step  toward  correcting  the  critical 
environmental  problems  facing  the  world  today.  Specifically,  one  of  the  problems 
addressed  was  chemical  products  damaging  the  protective  bubble  we  live  in.  This  bubble 
or  thin  layer  is  known  as  the  ozone  layer,  and  is  what  makes  the  earth's  atmosphere 
habitable  and  hence  unique  finm  any  other  known  planet  Without  diis  protective  layer,  the 
earth's  climate  would  not  support  life  as  we  know  it  today.  Protecting  this  ozone  layer 
from  chemical  destruction  is  paramount 

A  group  of  chemicals  that  deplete  the  ozone  layer  is  chlorofluorocarbons  (CFCs). 
CFCs  are  manmade  chemicals  of  chlorine,  flumine  and  carbon  and  were  developed  almost 
60  years  ago  by  General  Motors.  They  are  unique  in  that  they  have  a  combination  of 
desirable  properties  not  often  found  in  chemistry;  low  in  toxicity ,  non-flammable,  non- 
coTTOsive,  non-explosive,  extremely  stable  and  compatible  with  many  other  materials. 
However,  it  is  this  extreme  stability  which  causes  problems  to  the  ozone  layer,  because 
they  only  break  down  in  the  upper  atmosphere  (strato^heie)  when  subjected  to  intense 
ultraviolet  radiation.  This  breakdown  produces  chlorine  which  has  been  linked  to  the 
depletion  of  the  earth's  protective  ozone  layer.  A  chlorine  atom  destroys  an  ozone  molecule 
(O3)  by  reacting  with  it  and  converting  it  to  an  ordinary  oxygen  molecule  (O2)  and  carbon 
nxmoxide.  More  disturbing  is  the  fact  that  the  chlorine  atom  survives  the  ozone-destroying 
reaction  and  is  free  to  repeat  the  process  -  cxie  CFC  molecule  being  able  to  destroy  100,000 
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ozone  molecules.  Unfortunately,  a  typical  chlorine  molecule  can  survive  for  up  to  100 
years,  so  a  majcH*  problem  does  exist 

Armed  with  these  facts,  24  nations  (representing  the  United  Nations  Environment 
Program  (UNEP))  met  and  signed  the  Montreal  Protocol  in  1987.  They  discussed 
substances  that  deplete  the  ozone  layer  [Ref.  1]  and  called  for  a  near-term  freeze  on  the 
production  and  consumption  of  these  substances.  To  be  more  specific,  it  required 
production  of  these  chemicals  to  be  cut  back  to  1986  levels  followed  by  a  two-phased 
reduction  culminating  in  cutbacks  of  50%  by  mid-1998;  this  came  into  effect  on  July  1, 
1989.  In  1990,  a  progress  meeting  was  held  in  London  where  UNEP  delegates  agreed  to 
completely  phase  out  all  CFCs  by  the  year  2000  [Ref.  2].  This  was  a  needed,  but 
worrisome  change  for  the  U.S.  Navy  as  they  use  a  number  of  different  CFCs  for  various 
refrigeration  and  solvent/cleaning  needs.  One  such  need  is  the  approximately  1850 
shipboard  air  conditioning  (AC)  plants  presently  using  CFC-12  (in  reciprocating 
compressor  plants)  and  R-1 14  (in  centrifugal  compressor  plants).  These  AC  plants  not 
only  provide  personal  comfort,  but  cool  vital  electronic  con^)onent5,weapons  systems  and 
other  auxiliary  equipment 

To  comply  with  the  Montreal  Protocol  and  U.S.  legislation,  the  Mechanical  Systems 
Branch/(Code  2722)  at  David  Taylor  Research  Center  (DTRC)  is  pursuing  research  aimed 
at  eliminating  shipboard  use  of  CFC^s.  As  mentioned  by  Krinsky  and  Noel  [Ref.  3]  this 
research  is  to  be  completed  in  three  phases: 

1.  In  the  short  term,  to  identify  suitable  alternative  ozone-safe  chemicals  to  directly 
rqrlace  CFC-12  and  CFC-1 14.  In  order  to  do  this,  the  heat  transfn’  characteristics  must  be 
similar  to  the  existing  refrigerants  in  place  and  hence  the  need  for  a  database  exists  for 
current  refrigerants  (i.e.  R-1 14,  R-12,  R-1 13)  so  that  they  can  be  compared  to  the  new 
proposed  refrigerants. 


2.  In  the  longer  term,  to  research,  develop  and  test  substitute  chemicals  and 
alternative  technologies  to  replace  existing  CFC  uses. 

3.  To  implement  new  cooling  system  technologies  into  the  Fleet  which  do  not 
depend  on  CFCs  or  their  replacements.  One  such  technology  under  investigation  at  the 
Naval  Postgraduate  School  is  thermoacoustic  lehrigeration,  which  uses  sound  waves  to 
compress  and  expand  harmless  gas  molecules  (such  as  helium)  inside  a  tube. 

The  purpose  of  this  thesis  is  to  supplement  DTRCs  research  on  alternatives  to  CFCs 
by  establishing  baseline  nucleate  pool  boiling  data  of  R-113  from  a  small  bundle  of 
enhanced  tubes  representing  a  section  of  a  flooded  evaporator .  Emphasis  is  placed  on 
surface  history  affecting  activation/deactivation  of  nucleation  sites,  hysteresis,  and  the 
affect  of  varying  evaporator  pool  height  The  data  obtained  are  expected  to  save  as 
reference  data  for  comparison  with  future  proposed  refrigerants  from  bundles  of  both 
smooth  and  enhanced  tube  surfaces.  It  is  felt  that  what  is  learned  from  this  thesis  about  the 
charactoistics  of  R-1 13  can  be  applied  to  any  general  refrigerant  and  can  help  the  Navy 
design  future  AC  systems  using  alternative  chemicals. 

B.  OBJECTIVES 

The  objectives  of  this  thesis  are  as  follows: 

1.  Understand  in  greater  detail  both  the  convection  and  nucleate  pool  boiling 
phenomena  and  hysteresis  effects  within  a  small  enhanced  tube  bundle. 

2.  Vary  the  evfporator  pool  height  above  the  tube  bundle  to  see  how  this  effects 
nucleatimi  site  activation  and  general  evaporatm’  performance. 

3.  Investigate  the  effect  of  surface  history  on  nucleation  site  activation/deactivation. 

4.  Compare  data  with  similar  studies  (i.e.  same  fluid  or  tube)  done  at  the  Naval 
Postgraduate  School 
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II.  LITERATURE  SURVEY 


A.  GENERAL  INTRODUCTION 

Due  to  the  desirabili^  of  operating  many  engineering  devices  in  the  nucleate  boiling 
region,  many  types  of  enhanced  tubes  are  being  studied.  Nishikawa  et  al.  [Ref.  4]  noted 
that  there  are  two  primary  ways  to  increase  the  heat  transfer  performance  of  a  boiling  tube. 
The  first  was  to  decrease  the  wettability  of  the  fluid  in  which  you  are  immersed.  They 
suggested  using  a  Teflon  coating  on  the  boiling  surface.  However,  this  is  not  very 
effective  for  refrigerants  as  they  are  very  highly  wetting  m  any  surface  and  fill  the  cavities 
of  the  Teflon  surface.  The  second  way  proposed  was  to  manufacture  a  tube  with  many  re¬ 
entrant  cavities,  thereby  trapping  the  vapor  and  keeping  the  nucleation  sites  active.  These 
re-entrant  cavities  are  designed  such  that  if  one  site  is  activated,  neighboring  sites  will  soon 
follow  and  cause  the  entire  tube  to  nucleate  almost  instantaneously  due  to  the 
intercormecting  passageways.  There  are  now  many  such  tubes  with  re-entrant  cavities  on 
the  market 

The  two  main  types  of  reboilers  analyzed  today  are  kettle  reboilers  and  full  bundle 
boilers.  The  difference  is  described  in  detail  in  Payvar  [Ref.  S].  To  sununaiize,  in  the 
kettle  type  reboiler  the  tube  bundle  only  occupies  the  middle  lower  portion  of  the  evaporate 
leaving  large  areas  on  the  sides  of  the  bundle  for  fluid  circulation.  The  full  bundle  boilers 
occupy  the  whole  cross  section  of  the  lower  half  of  the  evaporator  shell.  As  can  be 
reasoned  from  these  definitions,  circulation  patterns  are  very  different  for  the  two  types  of 
boilers. 

The  following  is  a  brief  summary  of  the  work  done  on  the  bundle  apparatus  at  the 
Naval  Postgraduate  School: 
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1 .  TJ.  Murphy,  September  1987  [Ref.  Q 

•  Designed  System 

•  Conducted  Preliminary  Operadcm 

2.  CL.  Anderson,  June  1989  [Ref.  7] 

•  Tested  Smooth  Tubes,  R-113 

•  Tested  Smooth  Tubes,  R-1 14/Oil 

•  Tested  Finned  Tubes  (19  fyi),  R-1 14/Qil 

3.  N.  Akcasayar,  December  1989  [Ref.  8] 

•  Tested  Hnned  Tubes  (19  ^i),  R-1 14/Oil 

•  Tested  High  Flux  Tubes,  R-1 14/Oil 

4.  R  Eraydin,  December  1990  [Ref.  9] 

•  Tested  High  Flux  Tubes,  R-1 14/Oil 

•  Tested  Turbo-B  Tubes,  R-1 14/Oil 

B .  EXPERIMENTAL  AND  THEORETICAL  STUDIES 

In  the  literature  survey,  no  papers  wm  found  on  the  Turbo-B  tube  operating  in  a 
bundle  ot  how  surface  history  affects  Turbo-B  nucleation  site  activation/deactivaticni.  In 
addition,  no  information  was  found  on  varying  the  pool  height  above  the  bundle  to  see  how 
this  affects  nucleation  site  activation  and  general  evaporator  poformance.  These  areas 
could  be  very  inqxntant  to  the  U.S.  Navy  as  t)^ll  as  to  industry  and  could  cause  the  startup 
procedure  of  refrigeration  units  to  be  altered,  thus  emphasizing  the  need  for  the  objectives 
of  this  thesis. 

In  the  past  decade,  there  have  been  significant  improvements  made  to  tube  surfaces 
and  in  understanding  the  natural  convection  and  nucleate  boiling  regions  for  a  single  tube 
apparatus  and  several  correlations  have  been  adapted.  However,  the  tube  bundle  opens  a 
new  arena  of  topics. 
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Yilmaz  and  Westwater  [Ref.  10]  studied  the  effect  of  velociQr  past  a  smooth  tube  and 
how  it  effects  the  heat  transfer  performance.  They  used  a  6.5  mm  diameter  tube  in  R-1 13. 
They  found  that  with  an  increase  in  velociQr,  the  heat  transfer  performance  of  the  tube  was 
better. 

Cornwell  and  Schuller  [Ref.  11]  stated  that  the  increases  in  heat  transfer  coefficient 
found  by  Yilmaz  and  Westwater  could  not  oitirely  be  explained  by  the  increase  in  velocity 
past  the  tube.  Therefore,  they  decided  to  conduct  a  photographic  study  of  boiling  R>1 13  in 
a  bundle  at  one  atmosphere.  Their  study  used  241  electrically  heated  tubes,  19.1  mm  in 
diameter  and  arranged  in  a  square  in-line  pitch  oi  25.4  mm,  giving  a  pitch-to-diameter  ratio 
of  1.33.  They  found  that  the  bubbles  leaving  the  lower  tubes  in  the  bundle  impacted  and 
caused  a  sliding  motion  around  the  upper  tubes.  Cornwell  and  Schuller  attribute  the 
increase  in  the  upper  tubes  heat  transfer  performance  on  this  impact  and  sliding  motion  of 
bubbles  from  lower  tubes.  Cornwell  and  Schuller  also  showed  that  the  heat  transfer 
coefficient  was  the  highest  for  the  top  tube  in  the  bundle  and  decreased  as  one  nooved  down 
the  bundle;  the  difference  was  as  much  as  four  times.  Cornwell  [Ref.  12]  later  found  that 
in  the  nucleate  boiling  region,  sliding  bubbles  and  liquid  forced  convection  could  account 
for  all  the  heat  transfer  in  the  top  of  die  bundle. 

Chan  and  Shoukri  [Ref.  13]  studied  the  boiling  characteristics  of  a  small  enhanced 
multitube  bundle  in  a  pool  of  R-1 13.  The  rectangular  bundle  had  a  vertical  pitch  of  23.8 
mm  and  a  horizontal  pitch  of  31.75  mm.  From  there  studies,  they  concluded  that  the 
convective  contribution  is  highest  for  tubes  near  die  top  of  the  bundle  due  to  the  increased 
quality  and  fluid  velocity  across  these  tubes.  They  also  noted  at  lower  heat  fluxes,  the 
ccmvective  effects  dominate  whilst  for  higher  heat  fluxes,  the  primary  mode  of  heat  transfer 
was  nucleate  boiling.  At  high  heat  fluxes  thoe  are  numerous  active  nucleadon  sites  and  a 
high  bubble  generation  rate  which  decreased  the  influence  of  rising  bubbles  generated  from 
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tubes  below.  They  also  found  that  the  convective  effect  from  a  boiling  tube  below  an  upper 
tube  in  a  bundle  has  less  of  an  effect  as  it  gets  farther  from  the  tube  of  interest 

A  study  by  Bergles  et  al.  [Ref.  14]  showed  that  if  a  tube  is  on  the  verge  of  nucleation, 
it  burst  into  nucleation  if  tapped.  This  test  was  repeated  many  times  with  the  same  results. 
Several  other  experimenters  using  different  tubes,  refrigerants,  and  bundle  arrangements, 
have  also  shown  the  influence  of  tube  position  on  heat  transfer  performance  in  the  bundle. 
Wallner  [Ref.  IS]  used  a  smooth  tube  bundle  with  both  square  and  triangular  geometric 
arrangements  giving  pitch-to-diameter  ratios  of  1.33  to  1.5  respectively.  He  found  a  50 
percent  increase  in  the  heat  transfer  coefficient  over  a  single  tube  for  tubes  near  the  top  of 
the  bundle. 

Fujita  et  al.  [Ref.  16]  used  both  smooth  and  enhanced  tubes  to  study  the  effects  of 
tube  position,  distribution  of  heat  flux,  and  system  pressure  changes  (0.1  to  1  MPa)  within 
the  bundle.  They  used  a  triangular  tube  arrangement  with  R-113.  In  these  tests,  they 
reported  similar  results  as  above  with  higher  tubes  giving  improved  performance.  The 
enhancement  in  heat  transfer  again  was  believed  to  be  caused  by  bubbles  generated  from 
die  tubes  below. 

Webb  et  al.  [Ref.  17]  pointed  out  that  in  air  conditioning  applications,  the  refrigerant 
typically  enters  the  bundle  with  a  15  percent  quality.  With  this  entering  quality  into  the 
bundle,  it  simulates  the  effect  of  lower  tubes  activated  in  the  bundle  and  thus  will  inqirove 
the  performance  of  the  entire  bundle.  Furthomore,  the  boiling  mechanism  for  a  flooded 
evaporator  is  very  different  from  the  kettle  type  evaporator  and  therefore  it  is  very  hard  to 
formulate  and  evaluate  a  theoretical  model  for  a  bundle. 

Leong  and  Cornwell  [Ref.  18]  looked  at  the  same  representative  slice  of  a  flooded 
evaporator  that  Cornwell  and  Schuller  [Ref.  11]  and  Cornwell  [Ref.  12]  studied.  They 
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concluded  that  the  heat  transfer  coefficients  were  the  highest  at  the  top  of  the  bundle  and 
decreased  in  performance  as  you  went  down  the  bundle. 

Based  on  the  above-mentioned  studies,  all  experimenters  seem  in  agreement, 
independent  of  the  fluid  or  tube  used  in  the  bundle,  that  the  presence  of  boiling  lower  tubes 
causes  a  convective  flow  in  the  evaporator  that  gives  rise  to  an  increase  in  the  heat  transfer 
coefficient  cm  the  upper  tubes. 
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III.  EXPERIMENTAL  APPARATUS 


A.  TEST  APPARATUS  OVERVIEW 

The  whole  experimental  apparatus  including  the  auxiliary  equipment  and  the 
evaporator/condenser  test  apparatus  used  during  this  investigation  is  shown  in  Figure  1. 
The  following  is  only  a  general  description  of  the  whole  experimental  apparatus.  A  more 
detailed  locdc  at  the  condenser  and  evaporator  is  provided  in  Section  D.  Further  infnmadon 
about  the  apparatus  is  provided  by  Murphy  [Ref.  6]  and  Anderson  [Ref.  7]. 

Locddng  at  Figure  1,  the  apparatus  uses  three  closed  loops.  The  first  loop  ccxisists  of 
an  8  ton  re^geration  unit  located  outside  the  laborauvy  which  is  used  to  cool  an  ethylene 
glycol/water  mixture.  The  second  loop  is  the  ethylene  glycol/water  mixture  flowing 
through  the  condenser.  This  flow  rate  through  the  condenser  was  delivered  by  two  pumps 
which  could  be  operated  independently  ot  together,  this  coolant  mixture  condenses  the 
refrigerant  vapor  in  the  condenser  and  maintains  system  pressure.  Pump  number  one 
provided  coolant  flow  through  the  four  test  tubes  as  well  as  to  one  of  the  auxiliary  coils. 
Pump  number  two  provided  coolant  dirough  a  manifold  which  distributed  the  coolant  to  die 
remaining  four  coils  in  the  condenser.  The  third  loop  is  the  evaporator/condenser  itself. 
The  v^Kir  generated  in  the  evaporator  flowed  upward  and  condensed  in  the  condenser 
which  dien  returned  the  condensate  to  die  ev^xiraior  via  gravity. 

B.  AUXILURY  EQUIPMENT 

1 .  8  Ton  Refrigeration  Unit 

This  unit  was  used  to  cool  the  ethylene  glycol/water  mixture  (coolant)  to  the 
desired  temperature  needed  to  condense  the  refrigerant  vapor.  For  the  present  pool  boiling 
experiments,  the  temperature  control  was  set  to  maintain  the  sump  at  5**C.  This 
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refrigeration  unit  had  a  cooling  capacity  of  28  kW  (8  tons)  with  a  30  gpm  pump  providing 
circulation  of  the  coolant  to  achieve  the  set  temperature.  To  prevent  damage  of  the  unit, 
several  safety  features  were  provided:  low  temperature  cutout,  high  pressure  cutout,  low 
flow  cutout,  and  a  hot  gas  bypass  which  allowed  for  continuous  operation.  For  operation 
procedures,  see  Appendix  E. 

2.  Ethylene  GlycolAVater  Mixture 

The  coolant  used  was  a  54%  by  weight  ethylene  glycol/water  mixture.  This 
was  used  to  control  the  system  pressure  by  circulating  it  through  the  condenser. 

3 .  Pumps 

There  were  two  pumps  available  to  circulate  the  coolant  from  the  sump  through 
the  condenser.  Pump  number  one  fed  four  test  tubes  and  one  of  the  secondary  condenser 
coils.  Pump  number  two  was  used  to  feed  the  other  four  secondary  condenser  coils. 
Pump  number  one  provided  enough  cooling  for  tests  up  to  five  evaporator  tubes. 
Howeva,  during  bundle  and  simulation  evaporatcn-  operation,  both  pumps  were  needed  at 
heat  fluxes  above  approximately  60,000  W/rri2. 

4 .  Flowmeters 

Four  calibrated  float-type  flowmeters,  connected  to  pump  number  cme,  were 
used  to  measure  the  flowrate  passing  through  the  four  test  ctmdenser  tubes.  One  addititmal 
flowmeter  which  was  connected  to  pump  number  two  was  used  to  measure  the  total  flow  to 
the  secondary  condenser  coils.  In  addition,  each  of  the  four  secondary  condenser  coils  had 
a  globe  valve  to  further  regulate  (or  shut  off)  flow  as  desired.  In  the  test  tubes,  the  flow 
was  regulated  by  a  ball  valve. 

C.  DATA  ACQUISITION  SYSTEM/INSTRUMENTATION 

As  described  by  Akcasayar  [Ref.  8],  a  Hewlett  Packard  HP-3497A  Data  Acquisition 
System,  HP-9125  computer  and  HP-701  printer  were  used  for  data  acquisition,  data 
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reduction  and  data  printing  respectively.  Although  an  HP-9826  computer  and  HP-7470A 
plotter  can  be  used  for  final  graph  printing,  a  Macintosh  Classic  computer  using 
CricketGraph  1.3.2  was  utilized.  HP  Basic  3.01  was  used  for  data  reduction.  At  the 
beginning  of  every  set  of  runs,  the  HP-9125  computer  had  to  be  initialized  by  using  three 
system  discs  if  the  system  was  turned  off.  As  described  by  Anderson  [Ref.  7],  type-T 
copper-constantan  thermocouple  measurements  (mvolts)  were  made  on  the  HP-3497A  with 
the  relay  multiplexer  assembly  equipped  with  thermocouple  compensation.  A  20  channel 
relay  multiplexer  card  was  used  to  measure  the  voltage  output  from  voltage  and  amperage 
sensors.  Voltage  measurements  were  taken  from  separate  sensors  that  measured  the 
voltage  going  to  the  tube  bundle,  simulation  and  auxiliary  heaters.  The  total  amperage 
going  through  the  auxiliary  and  simulation  heaters  was  measured  using  an  American 
Aerospace  Control  (AAC)  current  sensor.  The  currents  of  each  instrumented  tube  heater 
were  measured  using  five  identical  current  sensors.  The  voltage  supplied  to  the  other  active 
tubes  was  also  measured  but  the  current  of  each  active  tube  was  not  (the  total  current  for  a 
pair  of  active  tubes  was  measured).  This  was  felt  to  be  sufficient  since  these  tubes  each  had 
the  same  power  output  (1000  W)  as  the  test  tube  heaters  and  there  was  no  apparent  reason 
to  rntmitor  each  active  tube  heat  flux  individually. 

Computer  channel  assignments  for  data  acquisition  and  array  assignments  are  given 
in  Table  1. 

D.  EVAPORATOR/CONDENSER 

An  overall  view  of  the  evaporator/condenser  is  shown  in  Hgure  2.  The  evaporator 
was  designed  to  simulate  a  small  portion  of  a  refrigerant  flooded  evaporator.  Front  and 
side  views  of  the  evaporator  are  shown  in  Figures  3  and  4.  It  was  fabricated  from  stainless 
steel  plate  and  formed  into  a  short  cylinder  610  mm  in  diameter  and  241  mm  long. 
Electrically  heated  tubes  were  cantilever-mounted  from  the  back  wall  of  the  evaporator  to 
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permit  viewing  along  the  axis  of  the  tubes  through  the  lower  of  two  glass/Plexiglass 
windows  mounted  in  the  front  The  glass  was  used  on  the  refrigo^t  side  in  order  to 
prevent  crazing  or  cracking  of  the  Plexiglass  and  the  Plexiglass  gave  the  glass  added 
strength.  Figure  5  is  a  schematic  sectional  view  of  the  evaporator  that  shows  the  four  kinds 
of  heated  tubes  used  in  this  study.  The  four  kinds  of  heater  tubes  were  as  follows: 
Instrumented,  Active,  Auxiliary,  and  Simulation.  Table  2  gives  the  power  rating  for  the 
heaters  and  the  number  used  in  the  evaporator. 

The  electric  power  can  be  applied  separately  to  each  set  of  heaters  using  a  STACO 
240V,  23.5  KVA  rheostat  controller  shown  in  Figure  6.  Also,  the  desired  number  of 
instrumented  tubes,  active  tubes,  auxiliary,  ot  simulation  heaters  can  be  turned  on  by  using 
individual  circuit  breakers.  The  four  auxiliary  heaters,  each  with  a  maximum  rating  of  4 
kW,  were  installed  two  on  each  side  of  the  test  bundle  to  maintain  the  liquid  pool  at 
saturated  conditions  and  to  provide  system  pressure  control.  The  five  simulation  heaters, 
also  capable  of  4  kW  each,  were  mounted  below  the  test  bundle  in  order  to  simulate 
additional  tube  tows  in  a  larger  bundle  and  to  provide  inlet  v^por  quality  into  the  bottom  of 
the  test  bundle  as  suggested  by  Webb  [Ref.  17].  The  auxiliary  and  simulation  heaters  had 
to  be  used  with  caution  to  avoid  breakdown  of  the  refiigerant 

The  tube  bundle  consists  of  instrumented,  active,  and  dummy  tubes.  The  location  of 
each  tube  is  represented  by  the  respective  letter  I,  A,  and  D  shown  in  Figure  S.  The  test 
bundle  consisted  of  two  types  of  heated  tubes;  active  tubes  (mariced  "A")  which  contained 
1  kW  cartridge  heaters,  and  instrumented  tubes  (marked'T')  which,  in  addition  to  the  1  kW 
cartridge  heaters,  contained  six  wall  thermocouples.  In  measuring  boiling  heat-transfer 
coefficients,  great  care  must  be  exercised  with  the  cartridge  heater  and  temperature 
measuring  instrumentation  to  ensure  good  accuracy.  For  example,  the  type  of  heater  used 
could  have  an  affect  on  heat  transfer.  Bergles  [Ref.  19]  staled  that  the  method  of  heating  is 
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iminateiial  except  for  the  case  of  a.c.  resistance  heating  which  may  cause  bubble  generadtm 
in  phase  with  the  power  supply.  However,  more  recently,  Jung  and  Bergles  [Ref.  20] 
who  carried  out  extensive  pool  boiling  data  with  R-113,  concluded  that  the  heat  transfer 
coefficient  of  a  single  tube  in  pool  boiling  is  not  sensitive  to  variations  in  the  cartridge 
heater  heat  flux  provided  that  enough  thermocouples  are  used  to  measure  an  average  wall 
temperature.  Prior  to  this  investigation,  the  instrumented  test  tubes  were  fabricated  in  a 
similar  way  to  those  used  by  Hahne  and  Muller  [Ref.  21]  and  Wanniarachchi  et  al.  [Ref. 
22].  The  exact  procedure  can  be  found  in  Eraydin  [Ref.  9].  Figure  7  is  a  cross-sectional 
sketch  of  an  instrumented  tube  used  during  this  thesis,  showing  the  construction  details  and 
the  location  of  the  wall  thermocouples.  The  thermocouples  were  imbeded  in  the  wall  at 
different  circumferential  and  longitudinal  positions  along  the  heated  section  of  the  tube. 
The  instrumented  tubes  were  located  along  the  centerline  of  the  tube  bundle,  forming  a 
vertical  in-line  column.  All  the  instrumented  and  active  tubes  were  Turbo-B  tubes  made  by 
Wolverine.  These  tubes  were  14.15  mm  in  diameter  and  were  arranged  in  an  equilateral 
triangular  pitch  (i.e.  centerline-to-centerline  spacing)  of  19.1  mm,  giving  a  pitch-to- 
diameter  ratio  of  1 .35.  The  bundle  also  contained  a  number  of  unheated  dummy  smooth 
tubes  (maiked  "D")  that  were  used  to  guide  the  two-phase  mixture  through  the  bundle.  The 
dummy  tubes  were  made  from  commercially  available  15.9  mm  in  diameter  smooth  copper 
tubing.  Two  vertical  baffle  plates  were  made  of  aluminum  to  restrict  circulation  into  and 
out  of  the  bundle  from  the  sides  .  A  small  open  space  approximately  5  mm  in  height  was 
left  between  the  bundle  and  the  dummy  tube  rack.  This  space  allowed  some  liquid  to  enter 
the  bundle  and  replace  the  vapor  being  generated  in  the  bundle.  However,  there  was  also  a 
space  below  the  dummy  rack  that  provided  the  majority  of  the  circulation.  Thus, 
liquid/vapor  circulation  was  vertically  upward  over  the  five  instrumented  test  tubes  with  no 
net  horizontal  component.  The  dummy  rack  (Figure  8)  used  12  tubes  made  of  solid 
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aluminum  rods  each  15.9  mm  in  diameter  and  spaced  19  mm  from  centoline  to  centerline. 
This  rack  had  a  triangular  pitch  arrangement  with  vapor  retaining  plates  on  the  side  and  was 
used  for  two  purposes.  This  dummy  rack  was  designed  to  collect  all  rising  two-phase  flow 
generated  by  the  simulation  heaters  and  direct  it  to  the  test  section.  The  other  purpose  was 
to  simulate  the  vapor  passing  though  a  larger  bundle  before  reaching  the  instrumented 
tubes.  The  liquid-vapor  mixture  after  passing  through  the  bundle  was  separated  when  it 
reached  the  pool  surface.  The  vapor  flowed  toward  the  condenser  and  the  liquid  was 
circulated  back  roward  the  bottom  of  the  bundle. 

The  condenser  included  four  test  tubes  (each  2  m  in  length)  in  a  vertical  in-line 
column  and  five  auxiliary  copper  coils  used  to  regulate  the  pressure  in  the  systenL  It  was 
designed  to  permit  independent  condensation  studies  of  small  in-line  tube  bundles,  using 
the  evaporator  as  a  source  of  vapor.  All  these  tubes  and  coils  were  cooled  by  a  refrigerated 
mixture  of  water  and  ethylene  glycol.  Each  of  the  test  tubes  and  auxiliary  coils  could  be 
operated  separately  using  ball  or  globe  valves  respectively  to  regulate  the  desired  flow  rate. 
The  vapor  that  was  produced  in  the  evaporator  was  guided  up  toward  the  top  of  the 
condenser  by  a  vapor  shroud  as  shown  in  Figure  9.  The  vapor  was  then  directed  down 
over  the  condenser  tubes  and  coils  where  it  condensed.  Once  condensed,  the  force  of 
gravity  returned  the  condensate  back  to  the  evaporator.  Observation  windows  in  both  the 
ccmdenser  and  the  evaporator  allowed  for  easy  viewing  of  condensation  and  boiling. 

E.  GEOMETRY  OF  ENHANCED  SURFACE  USED 

1.  Turbo-B  Tube 

The  Turbo-B  tube,  manufactured  by  Wolverine  Tube  Inc.,  has  a  enhanced 
surface  geometry.  The  exterior  boiling  enhancement  is  made  by  raising  integral  low  fins, 
cutting  diagonally  across  these  fins,  and  then  rolling  the  frns  to  compress  them  to  form 
mushroom-like  pedestals  [Ref.  23].  This  process  forms  numerous  re-entraint 
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passageways.  Figure  10  shows  the  surface  of  the  tube  at  25  times  its  actual  size.  The  tube 
is  currently  available  in  copper,  cupro-nickel,  and  low  carbon  steel.  The  relative 
dimensitms  of  the  tube  used  in  the  present  experiments  are: 

Tube  material  >  Copper 
Outside  diameter  =  14.15  mm 
Inside  diameter  » 12.7  mm 
Enhanced  surface  length  =  203.2  mm 
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Figure  1.  Schematic  View  of  the  Apparatus 
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Figure  2.  Evaporator/Condenser  Schematic 
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HIDDEN  VIEW  SHOWS  POSITION  OF  RACK  OF  DUMMY  TUBES 


178  mm  (  7!)00  PIPE 


Figure  4.  Side  View  of  Evaporator 
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FLANGE  STUDS  — - EAlmmOS") 

NOT  SHOWN 


Figure  5.  Sectional  View  of  Evaporator  Showing  Tube  Bundle 
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DUMMY  TUBE  RACK 


Figure  6.  Photograph  of  208  V,  75  A,  Variable  Transformers 
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Figure  7.  Thermocouple  Locations  on  an  Instrumented  Boiling 
Tube  and  Tube  Section  View. 
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Figure  8.  Photograph  of  Dummy  Tube  Rack 


Figure  9.  Vapor  flow  path  in  Condenser 


Figure  10.  Ciose-up  View  of  Turbo-B  Tube  Surface  (25  Times) 
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TABLE  1.  COMPUTER/DATA  ACQUISITION  ASSIGNMENT 


Channel 


00 


01 


02 


03 


04 


05 


40 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


51 


52 


T(10) 


T(11) 


T(12) 


Thermocouple  Description 


Vapor  1-Top  of  Condenser 


Vapor  2-Top  of  Condenser 


Vapor  3-Top  of  Evaporator 


Liquid  1-Top  of  bundle 


Liquid  2-Top  of  bundle 


Liquid  3-Bottom  of  bundle 


Tube  I.No.  1 


Tube  1,No.  2 


Tube  I.No.  3 


Tube  I.No.  4 


Tube  I.No.  5 


Tube  I.No.  6 


Tube  2.No.  1 


Tube  2.NO.  2 


Tube  2.No.  3 


Tube  2.No.  4 


Tube  2.NO.  5 


Tube  2.No.  6 


Tube  3.No.  1 


Tube  3.No.  2 
Tube  3.No.  3 
Tube  3.No.  4 
Tube  3.No.  5 
Tube  3.No.  6 
Tube  4.No.  1 
Tube  4, No.  2 
Tube  4.No.  3 
Tube  4.No.  4 
Tube  4.No.  5 
Tube  4.No.  6 
Tube  5.NO.  1 
Tube  5.No.  2 
Tube  5.No.  3 
Tube  5.No.  4 
Tube  5.No.  5 


I _ Tube  5,No.  6 _ I  69  |  T(35) 
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TABLE  1.  COMPUTER/DATA  ACQUISITION  ASSIGNMENT  (CONY.) 


Amperage  Sensor  Description 

Channei 

Array 

Tube  1 

30 

Amp(O) 

Tube  2 

31 

Tube  3 

32 

Tube  4 

33 

Tube  5 

34 

Amp(4)  I 

Active  Heater  Group  1 

35 

Active  Heater  Group  2 

36 

Amp(6)  I 

Active  Heater  Group  3 

37 

Active  Heater  Group  4 

38 

Active  Heater  Group  5 

39 

Auxiliary  Heaters 

25 

Simulation  Heaters 

26 

Voltage  Sensor  Description 

Channel 

Array 

Instrumented/Active 

27 

wmmsnmi 

Simulation  Heaters 

28 

Auxiliary  Heaters 

29 

TABLE  2.  EVAPORATOR  HEATERS 


Heater  Type 

Number 

Power  Rating  per  Heater 

Instrumented  Tube  Heaters 

5 

1000W 

Active  Tube  Heaters 

12 

1000W 

Auxiliary  Heaters 

4 

4000W 

Simulation  Heaters 

5 

4000W 

2‘ 


IV.  EXPERIMENTAL  PROCEDURES 


A.  REMOVAL  OF  TUBE  BUNDLE  FROM  EVAPORATOR 

Before  removing  the  tube  bundle  from  the  evaporator,  it  was  ensured  that  the 
evaporator  was  free  of  any  refrigerant  and  was  at  atmospheric  pressure.  The  tube  bundle 
support  block  bolts  were  then  loosened  and  removed.  The  tube  bundle  was  then  removed 
from  the  back  of  the  boiler.  Care  was  taken  when  removing  the  tube  bundle  as  to  not 
damage  the  thermocouple  or  heater  wires.  It  was  easier  to  disconnect  the  thermocouple 
wires  from  the  data  acquisition  board  and  the  heater  wires  from  the  power  panel  instead  of 
leaving  them  connected  to  the  tubes  when  it  was  necessary  to  work  on  or  clean  the  tubes  in 
the  bundle.  The  tube  bundle  consisted  of  five  instrumented  heated  Turbo-B  tubes,  10 
active  heated  Turbo-B  tubes,  18  dummy  smooth  tubes  and  two  dummy  Turbo-B  tubes. 
These  two  dummy  Turbo-B  tubes  were  equipped  with  heaters  and  could  be  conneaed  if 
desired  to  study  various  other  bundle  configurations  (see  Anderson  [Ref.  7]). 

B.  BUNDLE  DISASSEMBLY 

When  the  bundle  needed  to  be  disassembled,  it  was  ensured  that  there  was  a  clean 
woiidng  space  and  a  container  to  put  small  parts  in  during  the  procedure.  The  first  process 
was  to  remove  the  perspex  plate  attached  to  the  aluminum  baffle  plate  by  four  screws.  This 
plate  was  modified  such  that  the  eight  dummy  smooth  tubes  surrounding  the  instrumented 
tubes  extended  through  the  plate  (with  ptfe  plugs  which  held  the  tubes  in  the  proper  pitch- 
to-diameter  alignment).  This  alteration  became  necessary  since  at  high  heat  fluxes,  the 
smooth  tubes  (which  were  screwed  into  the  tube  support  block)  loosened  fiom  the  wall  and 
rattled  around.  The  ten  screws  on  each  aluminum  baffle  plate  were  then  removed.  The 
plates  were  then  pulled  off  the  bundle.  The  two  outside  smooth  durruny  tubes  remained  on 
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the  tube  bundle  support  block  as  they  were  countersunk  into  the  block.  The  three  inner 
smooth  tubes  could  be  removed  from  the  aluminum  block  as  they  were  attached  only  by  the 
screws  already  removed.  The  12  smooth  tubes  were  then  unscrewed  from  the  tube  bundle 
support  block  as  seen  in  Figure  1 1.  Care  was  taken  in  noting  the  order  in  which  all  the 
tubes  were  removed  as  it  saved  time  in  reassembly.  The  bundle  was  engraved  to  help  in 
this  process.  With  these  tubes  and  baffle  plates  removed,  only  the  instrumented  and  active 
heater  tubes  remained.  These  tubes  were  removed  by  loosening  the  outer  O-ring 
compression  plate;  this  allowed  the  tubes  to  be  pulled  from  the  tube  support  block.  Care 
was  taken  when  pulling  the  tube  from  the  tube  support  block  so  as  to  not  break  the 
thermocouple  or  heater  wires. 

C.  SYSTEM  CLEAN-UP 

If  the  system  was  contaminated  it  had  to  be  thoroughly  cleaned.  In  order  to  do  this, 
the  entire  apparatus  had  to  be  taken  completely  apart  and  cleaned  in  the  following  manner. 

After  removal  of  the  R-1 13  from  the  evaporator  and  system  at  atmospheric  pressure, 
all  electrical  connections  to  the  bundle  were  disconnected  and  the  front-viewing  glass 
windows  were  removed.  The  tube  bundle  was  then  be  removed  as  described  above.  The 
dummy  tube  rack  was  then  taken  out;  at  all  times,  the  system  was  ventilated  with  a  fan. 

Having  removed  the  tubes  from  the  tube  bundle,  they  were  individually  washed  with 
warm  water,  rinsed  and  then  wiped  down  with  acetone.  The  snKX>th  tubes  were  cleaned 
with  Copper  Brite  to  remove  any  tarnish.  They  were  wiped  down  with  warm  water  and 
acetone.  The  same  procedure  was  followed  for  the  Turbo-B  tubes  except  they  were  not 
cleaned  with  Copper  Brite.  During  the  cleaning  process,  a  soft  bristled  toothbrush  was 
used  to  ensure  the  enhanced  surface  was  cleaned  properly,  exercising  care  not  to  interfere 
with  the  tube  surface.  The  storage  tank,  if  thought  to  be  contaminated  also  was  drained, 
vented,  opened  and  cleaned  using  the  same  procedure. 
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D.  SYSTEM  REFIT 

1 .  Installation  of  tube  bundle  in  evaporator 

Prior  to  installing  the  evaporator  tubes  and  support  block,  the  bundle  and 
evaporator  was  thoroughly  wiped  down  with  acetone.  Once  the  evaporator  tubes  were 
installed  in  the  support  block,  and  prior  to  tightening  the  tube  bundle  support  block  nuts, 
the  bundle  was  guided  into  the  evaporator  secticai.  After  the  bundle  was  in  position,  it  was 
then  checked  that  the  dummy  rack  was  still  propo’ly  positioned  below  the  bundle  and  that 
the  vapor  thermocouple  positions  were  still  1.75  cm  above  the  bundle.  The  support  block 
nuts  were  then  tightened  equally  on  opposite  sides  to  give  equal  compression  of  the  gasket 
If  the  fix>nt  window  was  removed  during  maintenance,  the  window  was  put  back  on  very 
carefully  with  small  equal  torques  applied  to  each  nut  in  turn,  in  a  clockwise  direction 
around  the  outer  ring  support  After  the  window  was  in  place,  each  tube  (which  extended 
through  the  outer  O-ring  compression  plate)  was  tapped  forward  so  as  to  touch  the  firont* 
viewing-window.  Once  completed,  the  O-ring  compression  plate  was  tightened.  This 
compressed  the  O-iings  (which  were  slightly  undersized  because  the  tubes  were  so  close 
together)  providing  a  seal  between  the  tube  support  block  and  the  stainless-steel  backing 
plate.  Also,  the  compression  plate  had  grooves  for  the  O-rings  to  fit  into  to  ensure  proper 
alignment  and  seal. 

2.  System  Leakage  Check 

After  the  system  was  isolated  from  the  atmosphere  and  system  integrity 
restored,  a  Seargent  Welch  10  SCFM  vacuum  pump  was  turned  on  and  the  pressure  in  the 
system  was  reduced  to  25  in-Hg  vacuum.  The  system  was  then  secured  and  left  untouched 
for  at  least  10  hours  to  see  if  there  was  any  leakage.  If  there  was  significant  leakage(>l 
inHg),  then  the  system  vacuum  was  broken  by  cracking  valve  R-2  (Figure  2)  slowly  (to 
ensure  minimmn  moisture  was  allowed  to  entn*  the  system).  This  valve  was  left  open  until 
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the  system  was  at  atmospheric  conditions.  The  system  was  then  pressurized  with  air  to 
approximately  5  psig  through  valve  R-2.  Large  leaks  were  then  detected  by  simply 
listening  to  the  air  issuing  Grom  the  system.  If  the  leaks  were  small,  they  could  be  found  by 
applying  a  small  amount  of  soap/water  solution  to  all  joints  in  the  system.  One  leak  that 
was  found  was  where  the  sleeve  had  been  soldered  onto  the  Tuibo-B  tube.  Since  this  was 
found  to  be  the  source  of  a  leak,  all  the  Turbo-B  tube  joints  were  carefully  resoldered. 
After  all  leaks  had  been  corrected,  the  system  was  again  subjected  to  a  vacuum  for  a 
minimum  of  10  hours.  If  the  increase  was  less  than  1/2  to  1  inHg  rise  in  a  10  hour  period, 
then  the  system  was  filled  with  refrigerant 

E.  REFRIGERANT 
1.  Fill 

a.  From  System  Storage  Tank 

A  refrigerant  storage  tank  was  used  to  store  excess  R-113  during 
experiments.  The  storage  tank  prevented  discharge  of  the  R-113  into  the  atmosphere  and 
made  experimentation  less  costly.  To  fill  the  evaporator  with  R-1 13  from  die  storage  tank, 
the  system  pressure  was  equalized  to  the  stcnage  container  pressure  by  opening  valves  R-1 
and  R-7  (Figure  2).  Once  the  pressure  was  equal,  then  valves  R-6  and  R-4  were  opened. 
Due  to  the  pressure  difference  caused  by  the  hydrostatic  head  of  the  liquid  in  the  storage 
container,  R-113  flowed  from  the  storage  tank  into  the  evaporator.  The  amount  of 
refrigerant  that  was  transferred  was  controlled  by  throttling  valve  R-4  to  obtain  the  desired 
level.  Once  the  desired  amount  of  refrigerant  had  been  transfored,  valves  R-4,  R-6,  R-7, 
R-8  were  closed.  If  required,  additional  R-113  was  later  added  from  a  storage  drum  at 
atmo^heric  pressure  to  the  system  using  valve  R-2  (see  next  section). 
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b.  From  Refrigerant  Storage  Canister 

To  fill  the  apparatus  from  the  refrigerant  drum,  the  ethylene  glycol/watcr 
mixture  was  first  cooled  to  a  temperature  of  Both  condenser  pumps  were  put  into 
operation  thus  allowing  coolant  to  flow  through  the  condenser  section.  This  allowed  for  a 
differential  pressure  to  exist  between  the  evaporator  and  ambient  A  clean  plastic  hose  was 
connected  to  valve  R-2  (Figure  2)  and  the  other  end  was  end  put  in  the  refrigerant  drum. 
Once  this  differential  pressure  and  hose  hookup  was  completed,  refrigerant  transfer 
between  ambient  and  the  evaporator  was  possible  by  slowly  opening  valve  R-2. 

2 .  Removal 

a.  To  System  Storage  Tank 

For  tube  replacement  system  maintenance  or  clean-up  purposes,  the  R- 
113  was  transferred  to  the  storage  tank.  The  ethylene  glycol/water  mixture  temperature 
flowing  through  the  storage  tank  was  cooled  to  less  than  5^0,  valves  R-7  and  R-8  (Hgure 
2)  were  opened  and  the  vacuum  pump  was  turned  on  to  put  the  storage  tank  under  vacuum. 
Once  the  storage  tank  was  under  a  vacuum,  R-113  flowed  from  the  evaporator  to  the 
storage  tank  because  of  the  pressure  difference.  This  pressure  difference  was  maintained 
by  keeping  the  heaters  on  in  the  evaporator  and  the  coolant  flow  through  the  storage  tank. 
As  the  level  in  the  evaporator  came  down,  power  to  the  heaters  was  secured  to  avoid 
overheating.  Any  remaining  R-113  that  could  not  be  transferred  into  the  storage  tank  was 
drained  through  valve  R-5  and  put  into  a  drum  to  be  used  later. 

b.  To  Refrigerant  Storage  Canister 

This  procedure  is  primarily  used  to  change  refrigerants.  The  system 
vacuum  can  be  broken  by  slowly  cracking  open  R-2  (Figure  2)  to  ensure  that  moisture  was 
not  allowed  to  enter  the  system.  Once  the  system  was  at  atmospheric  pressure,  valve  R-2 
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was  opened  wide  to  allow  air  to  enter  when  the  refrigerant  was  removed.  The  refrigerant 
could  then  be  drained  through  valve  R-5  and  put  into  a  drum  to  be  used  later. 

F.  OPERATION 

1.  System  Lightoff,  Securing  and  Emergency  Procedures 

See  Appendix  E 

2.  Normal  Operation 

The  evaporatOT  was  filled  with  R-113  to  a  level  of  approximately  10  cm  above 
the  top  tubes  in  the  bundle.  Prior  to  operating  the  system,  the  8-ton  refrigeration  unit  was 
run  for  approximately  an  hour  to  reduce  the  ethylene-glycolAvater  mixture  in  the  sump  to  a 
temperature  of  S°C.  The  pressure  in  the  evaporator  was  usually  14  to  16  inHg  vacuum  if 
the  system  was  left  secured  overnight  As  the  sump  was  brought  to  temperature,  the  data 
acquisition  system  and  computer  were  turned  on.  This  allov/ed  the  temperature  in  the 
system  to  be  monitored  during  warm-up  to  saturation  conditions.  With  this  and  pump 
number  one  running,  the  auxiliary  heaters  were  then  used  to  bring  the  pool  up  to  saturatitm 
temperature.  Tlie  auxiliary  heaters  were  set  to  1  kW  nominal  after  pool  temperature  had 
been  reached.  The  system  pressure  was  slowly  raised  by  aUowing  minimal  coolant  to  flow 
at  a  low  rate  through  the  condenser  tubes  while  heating  the  system  with  the  auxiliary 
heaters.  Once  the  required  saturatitm  temperature  of  47.S*’C  was  reached,  die  instrumented 
tubes  were  switched  on  and  set  to  the  desired  value.  This  was  done  to  prevent  the  tubes 
from  nucleating  prematurely.  The  heat  flux  of  the  instrumented  tubes  was  then  slowly 
increased  by  adjusting  the  liieostat.  For  increasing  heat  flux,  the  data  were  taken  with 
very  small  heat  flux  increments,  waiting  at  least  S  minutes  to  attain  steady  state  conditions 
at  each  heat  flux.  At  all  regions  of  the  boiling  curve  (and  especially  near  the  onset  of 
nucleate  boiling),  two  readings  were  taken  for  each  heat  flux.  The  bundle  was 
ctmtinuously  monitored  through  the  observation  windows. 
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Figure  12  shows  the  tube  bundle  arrangements  used  during  experimentation. 
Test  One  was  with  only  instrumented  tube  number  one  on  in  the  bundle.  Test  two  was 
with  instrumented  tubes  one  and  two  on  in  the  bundle.  Test  three  was  with  instrumented 
tubes  one,  two  and  three  on  in  the  bundle.  Test  four  was  with  instrumented  tubes  one 
through  four  on  in  the  bundle.  Test  five  was  with  all  rive  instrumented  tubes  on  in  the 
bundle.  Test  six  was  all  five  instrumented  tubes  plus  all  five  pairs  of  active  heater  tubes  on 
in  the  bundle.  Test  seven  was  the  same  as  test  six  but  with  the  simulation  heaters  on. 

G.  DATA  REDUCTION  PROCEDURES 

The  data  reduction  program  "DRP4A''  was  used  during  the  experiments  for 
processing  the  data  collected.  The  program  was  written  in  HP  Basic  3.01  and  run  on  an 
HP-9000  series  computer.  The  characteristics  and  capabilities  of  this  software  and  the 
entire  listing  of  the  program  is  similar  to  that  provided  by  Anderson  [Ref.  7].  The 
following  modifications  were  made: 

1.  Conectitxi  for  pool  height 

2.  Installation  of  new  thermocouple  at  bottom  of  liquid  pool  (bundle  inlet 
ten:q)erature) 

The  entire  program  is  included  in  Appendix  F  to  provide  continuity  for  future  thesis 
work. 
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Figure  12.  Tube  Bundle  Arrangements  used  During  Experimentation 
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V.  RESULTS  AND  DISCUSSION 


A.  GENERAL  COMMENTS 

The  results  are  presented  in  five  sections  with  sub-sections  as  appropriate.  The  first 
section  discusses  the  natural  convection,  nucleate  pool  boiling  phenomena  and  hysteresis 
effects  within  a  small  enhanced  tube  bundle.  The  second  section  discusses  how  evaporator 
pool  height  above  the  tube  bundle  affects  the  nucleation  site  activation  and  general 
evaporator  performance.  The  third  section  looks  at  the  effects  of  auxiliary  heaters.  The 
fourth  section  discusses  the  effect  of  surface  history  on  nucleation  site 
activation/deactivation.  The  final  section  shows  comparisons  of  data  taken  during  this 
thesis  with  previously  obtained  data  at  the  Naval  Postgraduate  School 

A  list  of  data  runs  conducted  during  this  investigation  may  be  found  in  Appendix  A. 
AH  data  files  used  in  this  thesis  use  the  following  filename  sequence.  Each  file  is 
composed  of  five  sets  of  alpha-numeric  characters  used  to  describe  the  experiment 

First  set  (1  char.)  Heat  flux  I(Increasing  heat  flux) 

Dfl^ecreasing  heat  flux) 

Second  set  (2  char.)  Type  of  tube  used  TB  (Turbo-B) 

Third  set  (2  char.)  Number  of  Experiment 

for  accounting  purposes 

To  give  an  example,  the  file  name  rrB12  means  "increasing  heat  flux,  using  a  Turbo- 
B  tube  and  the  12th  data  run.  If  trxne  detail  is  desired  about  a  data  run,  see  Appendix  A. 
All  plot  files  are  similar  to  data  files  except  drey  start  with  the  letter  "P".  All  tests  were 
dtme  with  the  auxiliary  heaters  at  either  0, 1,  or  3  kW  nominal  power  setting  (for  details  see 
Appendix  A). 
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All  the  graphs  show  heat  flux  (W/m^  along  the  ordinate  (y  axis)  and  wall  superiieat 
(K)  along  the  abscissa  (x  axis).  The  heat  flux  was  changed  through  the  use  of  a  variac 
firom  600  to  100,000  W/m^  for  increasing  or  100,000  to  600  W/m^  for  decreasing  runs. 
Approximately  25-30  points  were  taken  on  the  increasing  heat  flux  and  about  20  points 
were  taken  for  the  decreasing  heat  flux.  The  wall  superheat  is  the  difference  between  the 
corrected  average  tube  wall  temperature  and  the  liquid  saturadtm  temperature. 

To  clarify  the  operation  of  the  evaporator  during  the  above  tests,  we  used  seven  basic 
bundle  arrangements  during  the  test  procedures  as  discussed  and  shown  in  Chapter  IV. 

Anderson  [Ref.  7]  used  four  types  of  surflu;e  preparation  techniques. 

1.  Surface  preparation  A  maintains  the  evaporator  tube(s)  at  a  heat  flux  of  30  kWAn^ 
for  one  hour,  slowly  reducing  heat  flux  to  a  minimum  value  of  around  1  kW/m^,  and 
beginning  data  collection  immediately  at  successively  increasing  heat  fluxes.  The 
motivation  for  this  surface  preparation  is  a  widely  varying  dynamic  load  often  encountered 
in  an  air-condidcxiing  machine. 

2.  Surface  preparation  B  simulates  an  operating  refrigeration  system  undergoing 
frequent  off-cycles.  A  heat  flux  of  30  kW/n^  at  saturation  conditions  is  maintained  for  one 
hour.  Power  is  then  secured  for  30  minutes  and  then  data  are  taken  at  successively 
increasing  heat  flux  steps. 

3.  Surface  preparation  C  represents  the  initial  start-up  of  a  refrigeration  system. 
During  this  process,  the  evaporator  was  brought  up  to  temperature  and  pressure  with  the 
auxiliary  heaters  alone  to  ensure  the  tubes  did  not  have  any  active  nucleation  sites.  Once  at 
saturation  pressure,  the  tube  was  turned  on  and  the  heat  flux  varied  from  600  to  100,000 
W/m^  to  develop  a  heat  flux  versus  wall  superheat  plot  This  process  is  very  interesting  as 
it  introduces  many  ctxxiplexities  such  as  boiling  hysteresis. 
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4.  Surface  preparation  D  represents  a  continuous  operating  air-conditioning  system. 
It  is  closely  modelled  by  a  continuously  decreasing  heat  flux  starting  from  maximum  power 
of  100  kW/m2  after  it  has  been  boiled  at  this  power  for  one-half  hour.  Considering  that  the 
return  flow  to  the  boiler  is  normally  two-phase,  nucleation  sites  are  maintained  active 
throughout  the  evaporator  tube  bundle.  The  decreasing  heat  flux  data  runs  are  ncmnally  the 
most  readily  reproducible  and  are  used  most  often  in  measurements  of  tube-bundle 
performance.  This  surface  preparation  is  the  most  widely  used  for  reporting  the 
peifomiance  of  evaporator-tube  surfaces. 

Three  of  these  surface  preparation  techniques  were  used  during  the  course  of  the 
present  wotic.  Surface  preparation  technique  C  and  D  were  used  exactly  as  stated  above 
and  surface  preparation  technique  B  was  the  same  but  with  a  moditication.  Surface 
preparation  technique  B  was  modiHed  such  that  power  to  the  tube  heaters  could  be  secured 
for  any  amount  of  time  but  the  pool  temperature  had  to  be  maintained  at  saturation 
temperature. 

1 .  Repeatability  of  data 

Figure  13  presents  similar  data  taken  on  two  different  days  for  the  top  tube 
operating  alone  in  the  bundle;  the  repeatability  is  seen  to  be  very  good.  The  arrows  on  the 
figure  indicate  the  familiar  hysteresis  loop.  The  two  tests  show  very  similar  temperature 
overshoot  values  (indicated  by  the  horizontal  line).  The  point  of  incipience  (point  where 
nucleate  boiling  first  occurs)  is  also  very  similar  fOT  both  tests. 

The  heat-transfer  coefficient,  h,  is  often  used  as  a  measure  of  tube  performance: 
it  is  found  by  dividing  the  heat  flux  by  the  wall  superheat  From  Figure  13,  values  of  h  in 
the  nucleate  boiling  region  range  from  1540  W/m^K  at  low  heat  flux  to  13,750  WAn^K  at 
high  heat  flux.  Corresponding  values  of  h  in  the  convection  region  range  from  480 
W/m^K  at  low  heat  flux  to  700  W/m^K  at  medium  heat  flux.  It  can  be  seen  from  Figure 
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13  that  the  tube  performance  typically  improves  by  a  factor  of  five  once  the  tube  starts  to 
nucleate.  This  explains  the  importance  of  ensuring  a  flooded  evaporator  operates  in  the 
nucleate  boiling  rather  than  the  natural  convection  region. 

Also  shown  on  Figure  13  are  typical  uncertainty  bands  for  the  data  at  both  high 
and  low  heat  flux  and  wall  superheat  It  can  be  seen  that  the  greatest  uncertainty  occurs  at 
low  wall  superheat  and  heat  flux;  this  was  as  expected  since  the  values  of  wall  superheat 
are  of  the  order  accuracy  for  the  thermocouples. 

Figure  14  shows  the  limited  predictions  that  are  available  for  both  natural 
convection  and  nucleate  pool  boiling  data.  These  correlations  were  formulated  for  single 
smooth  tubes,  but  represents  all  that  is  available  for  enhanced  surfaces  and  tube  bundles. 
This  indicates  the  urgent  need  for  further  work  in  this  area.  For  the  natural  convection 
region,  two  correlations  are  shown:  Churchill  and  Chu  [Ref.  24]  and  Churchill  and  Usagi 
[Ref.  25] .  The  experimental  data  for  a  single  tube  within  the  bundle  show  an  appreciable 
higher  heat-transfer  perfonnance  than  prediction.  This  might  be  expected  as  the  tube  is  an 
enhanced  surface.  However,  in  the  convection  region,  one  should  not  expect  a  large 
difference  in  performance  between  different  tube  surfaces  as  they  are  designed  primarily  for 
the  nucleate  region.  Indeed,  experiments  by  Murphy  [Ref.  6]  have  shown  that  the  data  for 
enhanced  surfaces  falls  on  top  of  smooth  tube  data  in  the  convection  region. 
Consequently,  die  discrepancy  seen  is  more  likely  to  be  due  to  circulation  effects  within  the 
bundle. 

In  the  nucleate  boiling  region,  there  is  also  a  large  enhancement  seen  in  the 
experimental  data  when  compared  with  the  correlation  of  Stephan  and  Abdelsalam  [Ref. 
26].  However,  in  this  region,  the  discrepancy  can  be  explained  due  to  the  enhanced  nature 
of  the  surface.  Circulation  effects  are  still  likely  to  be  present,  but  small  for  a  single  tube 
when  compared  to  the  large  number  of  active  nucleation  sites. 
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B .  BUNDLE  TRENDS  FOR  A  POOL  HEIGHT  OF  10  CM 

Rgure  15  shows  the  influence  of  tube  two  on  tube  one.  One  may  expect  the  thennal 
currents  from  the  lower  tube  to  further  enhance  the  upper  tube.  However,  with  this  lower 
tube  heated,  the  performance  of  tube  one  in  the  convection  region  did  not  change.  More 
striking  is  the  significantly  poorer  performance  of  tube  two  in  this  region  for  which  no 
explanation  can  be  found.  It  is  interesting  to  note  that  tube  two  (which  should  behave  as 
an  isolated  single  tube  in  the  convection  region,  i.e.  feeling  no  effect  from  the  tube  above) 
is  now  in  much  better  agreement  with  the  natural  convection  correlations. 

Also  seen  is  that  the  top  tube  nucleates  at  a  lower  heat  flux.  However,  once  both 
tubes  are  nucleating,  the  agreement  is  very  close,  with  the  top  tube  having  the  slightly 
better  performance.  Although  this  slight  discrepancy  could  be  due  to  uncertainty  in  the 
data,  it  is  more  likely  to  be  due  to  the  influence  of  the  bubbles  coming  from  the  tube  below, 
especially  when  one  con^aies  it  with  the  top  tube  acting  alone  (Figure  13). 

Figure  16  shows  decreasing  heat  flux  for  two  tubes.  Agreement  is  very  close,  with 
the  top  tube  performing  slightly  better  at  high  heat  flux,  but  the  opposite  at  low  heat  flux. 
However,  as  heat  flux  is  reduced,  the  uncertain^  in  die  data  increases. 

With  three  tubes  in  operation  (Figure  17) ,  the  trends  are  much  the  same.  One  point 
of  interest  is  the  effect  on  lower  tubes  as  upper  tubes  start  to  nucleate.  As  tube  one  starts 
to  nucleate,  tubes  two  and  three  shift  significantly  to  the  left,  indicating  an  increase  in 
heat-transfer  coefficient  However,  tubes  two  and  three  remain  in  the  convection  region. 
Any  further  noovement  to  the  left  by  the  upper  tubes  causes  a  corresponding  shift  to  the  left 
for  the  tubes  below.  This  is  probably  due  u>  increased  circulation  effects  (i.e.  as  the 
bubbles  from  the  upper  tubes  move  around  the  evaporator,  they  help  to  agitate  the  thermal 
boundary  layer  around  the  tubes  below).  The  shift  to  nucleation  is  strictly  in  order  as  one 
nooves  down  the  bundle. 
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For  decreasing  heat  flux  (Figure  18),  there  is  a  significant  drop-off  in  performance 
for  tube  three,  especially  at  a  low  heat  flux.  It  may  be  that  as  one  moves  down  the  bundle, 
the  greater  pressure  deactivates  the  nucleadon  sites  at  a  faster  rate. 

Figures  19  to  22  show  four  and  five  tubes  activated  for  increasing  and  decreasing 
heat  flux.  Similar  trends  to  those  already  reported  are  seen;  tube  one  nucleates  first 
followed  by  two,  three,  etc.  As  each  successive  tube  nucleates,  the  performance  of  the 
remaining  lower  tubes  is  improved  slightly.  In  the  nucleate  region  at  high  heat  flux,  tube 
one  exhibits  the  best  performance  followed  by  tubes  two,  four,  five  and  three.  This 
somewhat  strange  order  maybe  due  to  uncertainty  at  the  fact  that  tube  three  is  only  partially 
nucleating.  At  the  very  highest  heat  flux,  all  five  tubes  seem  to  coincide.  At  low  heat 
fluxes,  tube  two  gives  the  best  heat  transfer;  it  is  unclear  as  to  why  this  should  be. 

Because  of  the  unexpected  behavior  of  tube  three  at  high  heat  flux,  the  positions  of 
tubes  three  and  five  were  switched  to  see  if  this  behavior  was  a  "bundle"  effect  or  simply 
due  to  the  particular  tube  itself.  Figure  23  shows  the  results  after  switching  for  decreasing 
heat  flux.  New  tube  five  (fnmerly  tube  duee)  now  has  the  poorest  performance,  indicating 
that  it  is  the  tube  rather  than  the  bundle  causing  the  eariier  effect  It  may  be  that  air  pockets 
surround  the  thermocouple  junctions  of  tube  five.  Consequently,  the  tubes  were  kept  in 
this  positirm  fw  the  remainder  of  the  tests,  with  dte  "bad"  tube  in  the  lowest  position. 

Figures  24  and  25  show  the  entire  bundle  activated  for  both  increasing  and  decreasing 
heat  flux  while  Hgures  26  and  27  show  die  entire  bundle  plus  simulation  heaters  activated 
for  increasing  and  decreasing  heat  flux.  These  simulation  heaters  give  the  bottom  of  the 
bundle  an  inlet  quality  of  about  20  percent,  thereby  simulating  a  real  evaporator.  Similar 
trends  to  those  seen  before  are  evident  The  heat  flux  at  which  incipience  occurs  continues 
to  decrease  for  all  tubes,  indicating  the  extra  quality  moving  through  the  bundle  helps  to 
trigger  "early"  nucleation.  As  this  simulates  pan  of  a  real  evaporator,  this  trend  is 
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encouraging,  as  it  is  more  likely  to  be  operating  in  the  nucleate  region  rather  than  the 
convection  region.  For  the  decreasing  heat  flux  run,  the  heat  transfer  from  the  lower  tubes 
has  also  been  improved  due  to  the  increased  vapor  passing  over  them. 

As  already  mentioned,  it  is  difficult  to  make  too  many  direct  comparisons  when  all 
five  tubes  are  independently  instrumented.  One  of  the  tubes  may  have  certain  induced 
errors  due  to  manufacture  that  another  tube  does  not  Consequently,  of  perhaps  more  use 
is  to  look  at  what  happens  to  the  top  tube  only  as  successive  tubes  are  activated  below. 
Figure  28  does  this  for  all  seven  tests  (see  Chapter  IV:  nmmal  operation,  fen-  description  of 
tests)  for  an  increasing  heat  flux.  In  the  natural  convection  region,  the  point  of  incipience 
occurs  at  increasingly  lower  heat  fluxes  as  nrore  tubes  are  activated.  In  the  nucleate  boiling 
region,  the  heat-transfer  coefficient  increases  as  more  tubes  are  activated.  Both  the  natural 
convection  and  nucleate  boiling  regions  are  well  defined  indicating  good  reliability  of  the 
data.  For  decreasing  heat  flux  (Figure  29),  the  heat-transfer  coefficient  is  again  clearly 
increased  as  more  tubes  are  activated  below.  This  is  more  evident  at  low  heat  flux  than  at 
high  heat  flux. 

C.  BUNDLE  PERFORMANCE  WITH  VARYING  POOL  HEIGHT 

One  of  the  objectives  of  this  thesis  was  to  see  the  effect  of  pool  height  on  heat- 
transfer  performance.  Pool  height  here  is  defined  as  the  height  from  the  free  liquid  surface 
down  to  the  top  of  the  bundle  (tube  one).  All  experiments  to  date  had  been  carried  out  at  a 
pool  height  of  10  cm.  Two  additional  heights  of  0  cm  and  20  cm  were  also  studied  with 
the  same  vapor  pressure  within  the  evaporator  (i.e.  to  see  the  effect  of  pressure  head  on  the 
tubes  in  the  bundle).  It  is  understood  that  a  true  flooded  evaporator  has  zoo  pool  height 
(indeed  the  top  tubes  operate  in  the  bubbly  froth). 

Figures  30  to  42  show  the  same  seven  tests  for  a  pool  height  of  20  cm  as  previously 
shown  for  10  cm.  The  trends  of  all  these  graphs  are  exactly  the  same  as  for  10  cm.  The 
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effect  of  the  greater  pool  height  is  generally  to  delay  the  point  of  incipience.  Figures  43  and 
44  show  tube  one  only  for  all  seven  tests  (i.e.  to  see  the  effect  of  successive  activated  tubes 
below)  for  increasing  and  decreasing  heat  flux  respectively.  This  to  shows  similar  trends 
to  before,  with  the  effect  of  successive  tubes  triggaing  earlier  nucleation  on  the  top  tube  as 
well  as  enhaixnng  die  nucleate  boiling  heat-transfer  coeffident 

Figures  45  to  57  show  the  same  seven  tests  for  a  0  cm  pool  height  and  again  the 
trends  are  very  similar.  However,  this  time,  the  point  of  incipience  occurs  at  the  lowest 
value  of  heat  flux.  Figures  58  and  59  show  tube  one  only  for  all  seven  tests  for  increasing 
and  decreasing  heat  flux  respectively.  The  conclusions  drawn  from  these  figures  are  the 
same  as  before. 

To  summarize  the  effect  of  pool  height.  Figures  60  and  61  show  tube  one  (i.e.  Test 
One  )  at  the  three  pool  heights  for  increasing  and  decreasing  heat  flux.  It  can  be  seen  that 
the  zero  pool  height  data  nucleates  first  followed  by  10  cm  and  then  20  cm.  Similar  trends 
were  seen  regardless  of  what  test  was  chosen  i.e.  for  test  seven  with  whole  bundle  plus 
simulation  heaters  activated,  all  duee  points  of  indpience  (for  die  3  pool  hdghts)  move  to  a 
lower  heat  flux  but  the  zero  pool  height  data  still  nucleate  first  Obviously,  the  greater 
pressure  heads  felt  by  the  tubes  at  larger  pool  heights  inhibit  the  transition  to  nucleate 
boiling  i.e.  largo*  wall  superheats  are  required  to  overcome  the  increased  pressure.  Further 
experiments  should  be  conducted  in  the  future  whereby  the  vapor  pressure  above  the 
bundle  is  varied  with  the  pool  height  such  diat  the  pressure  diat  each  tube  "feels"  at  a  given 
locatitm  is  the  same. 

D.  EFFECT  OF  AUXILIARY  HEATERS 

For  all  tests  reported  so  far,  the  auxiliary  heaters  were  set  to  IkW  nominal.  The 
auxiliary  heaters  were  needed  to  maintain  pool  temperature  throughout  the  test  They  were 
only  really  needed  at  low  heat  fluxes  or  if  only  a  few  tubes  in  the  bundle  were  activated. 
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For  these  tests,  if  they  were  not  used,  the  evaporator  could  not  produce  enough  vaptnr  to 
keep  the  pressure  constant  even  with  the  condenser  shut  down  (i.e.  the  walls  of  the  test 
section  acted  as  the  condenser).  As  a  consequence,  the  pool  temperature  would  drop. 
However,  the  necessary  operation  of  the  auxiliary  heaters  has  an  effect  on  the  circulation 
patterns  within  the  bundle  (this  has  been  demonstrated  in  the  single  tube  facility). 
Therefore,  by  having  these  auxiliary  heaters  set  to  the  same  value  for  each  test,  any  effect 
they  may  have  should  be  a  constant  factor.  To  see  the  effect  that  the  auxiliary  heaters  have. 
Figure  62  shows  data  for  three  tubes  on  during  a  decreasing  heat  flux  with  the  auxiliary 
heaters  off  compared  to  Figure  18  with  them  %t  at  1  kW  for  a  10  cm  pool  height. 
Essentially  the  effect  is  minimal  in  the  nucleate  region  at  high  heat  fluxes.  At  lower  heat 
fluxes,  no  data  could  be  taken  with  the  auxiliary  heaters  off.  In  the  convection  region 
(again  at  low  heat  flux  so  there  is  no  data  with  the  auxiliary  heaters  ofO,  the  effects  of 
increasing  the  auxiliary  heater  power  is  to  increase  the  heat-transfer  .coefficient  (i.e. 
reducing  the  wall  superheat).  The  same  effea  was  seen  for  four  tubes  with  the  auxiliary 
heaters  off  and  set  to  1  kW  for  a  10  cm  pool  height;  shown  in  Figures  63  and  20 
respectively.  Figures  64  and  30  show  the  same  trend  for  tube  one  with  the  auxiliary 
heaters  set  at  3  kW  compared  to  1  kW  for  a  pool  height  of  20  cm.  These  trends  are 
probably  seen  due  to  the  increased  circulation  within  the  bundle  at  the  higher  auxiliary 
heatCT  setting.  Obviously,  it  would  be  preferable  to  do  away  with  the  auxiliary  heaters 
altogether.  With  R-1 14,  due  to  the  lower  normal  boiling  ptnnt,  this  should  be  possible. 

E.  NUCLEATION  SITE  ACTIVATION/DEACmVATION 

Nucleation  site  activation/deactivation  is  not  well  understood.  Once  a  tube  is 
nucleating,  there  are  a  coiq)le  of  ways  that  the  nucleation  sites  can  die  out: 

1.  Reducing  the  heat  to  the  tube  surface  thereby  allowing  the  fluid  to  collapse  the 
vapOT  bubble  and  re-wet  the  surface  due  to  pressure  effects. 
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2.  Increasing  the  subcooling  of  the  liquid  around  the  tube  (i.e.  causing  the  vapor 
bubble  to  condense  into  the  fluid). 

The  following  tests  were  conducted  to  see  how  quickly  the  nucleation  sites,  once 
activated,  died  away  if  the  power  to  the  tube  was  cut  off  (i.e.  to  see  how  quickly  the  liquid 
re-wets  the  surface  as  discussed  in  mechanism  one  above).  The  pool  temperature  was 
maintained  at  saturation  conditions  throughout  the  test  (i.e.  there  was  no  subcooling).  All 
the  tests  were  conducted  on  the  top  tube  only.  The  test  procedure  used  was  the  modified 
surface  preparation  B.  The  power  to  the  heater  was  then  turned  off  for  a  known  finite  time 
and  then  reactivated.  An  increasing  run  was  then  carried  out  This  was  repeated  fcv  a 
number  of  different  times  and  pool  heights.  Hie  auxiliary  power  setting  was  1  kW  for  all 
tests. 

Figures  65  and  66  show  the  data  for  a  pool  height  of  20  cut  The  turn-off  time  for  die 
heater  varied  between  five  seconds  and  one  hour.  It  can  be  seen  that  after  only  15  s,  the 
nucleation  sites  have  all  died  out  and  the  tube  behaves  as  if  it  had  been  secured  overnight 
(these  are  the  data  sets  listed  as  14  and  15  hours).  Of  course,  by  securing  overnight, 
subcooling  also  plays  a  part  in  deactivating  the  sites. 

Figure  67  shows  the  data  for  a  pool  height  of  10  cm.  The  turn-off  time  for  the  heater 
varied  between  five  seconds  and  one  minute.  It  can  be  seen  that  after  one  minute,  the 
nucleation  sites  have  all  died  out  and  the  tube  behaves  as  if  it  had  been  secured  overnight 
Hence,  it  would  seem  that  due  to  the  reduced  pressure  head  trying  to  push  the  liquid  into 
the  sites,  tiie  time  needed  to  deactivate  all  the  sites  is  significantly  Iraiger. 

Figures  68  and  69  show  the  data  for  a  pool  height  of  0  cm.  The  turn-off  time  for  the 
heater  varied  between  40  seconds  and  8  hours  (keeping  pool  temperature  constant 
throughout).  It  can  be  seen  that  even  for  a  shut-off  time  of  8  hours,  none  of  the  nucleation 
sites  are  deactivated.  Hie  pressure  head  is  negligible  here  and  it  is  mechanism  two  that 
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deactivates  the  sites.  As  the  pool  temperature  was  not  subcooled  during  this  particular  test, 
none  of  the  sites  are  therefore  deactivated.  Hgure  70  shows  a  test  where  the  apparatus  was 
secured  for  2.75  hours  (i.e.  allowed  to  subcool).  The  deactivation  seen  is  due  to  the 
subcooling  that  has  occurred. 

Similar  trends  for  all  three  pool  heights  have  been  obtained  for  an  auxiliary  heater 
power  setting  of  3  kW.  Slight  differences  in  the  times  needed  for  total  deactivation  were 
noted  to  be  20  seconds  for  10  cm  vice  one  minute.  These  data  are  shown  in  Figures  71  to 
75. 


F.  COMPARISONS  WITH  PREVIOUS  NPS  DATA 

Figure  76  shows  a  comparison  between  tiie  present  data  for  a  Turbo-B  bundle  in  R- 
113  and  the  data  of  Anderson  [Ref.  7]  for  a  smooth  tube  bundle  in  R-113.  For  clarity, 
only  the  top  tube  has  been  shown  for  a  decreasing  heat  flux.  It  can  be  seen  that  the  heat 
transfer  enhancement  given  by  the  enhanced  tube  is  three  to  five  times  for  high  and  low 
heat  fluxes  respectively. 

This  is  also  true  for  the  overall  bundle  heat  transfer  coefficient,  giving  a  factor  of  4 
enhancement  at  a  typical  operating  heat  flux  of 20,000  W/m^.  Also  shown  on  the  figure  is 
the  correlaticm  of  Stephan  and  Abdelsalam  [Ref.  26]  showing  a  fairly  good  agreement  witii 
the  smooth  tube  data. 

Figure  77  shows  a  comparison  between  the  present  data  and  that  of  Eraydin  [Ref.  9] 
who  used  a  Turbo-B  bundle  in  R-114.  It  can  be  seen  that  there  is  very  little  difference 
between  the  two  sets  of  data.  The  R-1 14  data  giving  the  slightly  better  performance.  This 
could  be  due  to  the  slightly  better  thermal  conductivity  of  R-114  at  its  saturation 
teiiq)erature  con^ared  to  R-1 13  at  its  saturation  temperature. 
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Figure  13.  Performance  of  Tube  1  at  a  10  cm  Pool  Height  for 
Similar  Data  Taken  on  Different  Days  to  Show  Repeatability 
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Heat  Flux/(W/m*2) 


Figure  14.  Performance  of  Tube  1  at  a  10  cm  Pool  Height  Along  With 
Correlations  for  Natural  Convection  and  Nucleate  Pool  Boiling  Regions 
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Figure  19.  Performance  Variation  of  Tubes  1  to  4  at  a  10  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB07) 
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Figure  20.  Performance  Variation  of  Tubes  1  to  4  at  a  10  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB07) 
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Figure  22.  Performance  Variation  of  Tubes  1  to  5  at  a  10  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB27) 
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Figure  25.  Performance  Variation  of  the  Bundle  at  a  10  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB48) 
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Figure  26.  Performance  Variation  of  the  Bundle  plus  Simulation 
Heaters  at  a  10  cm  Pool  Height  During  a  Increasing  Heat  Flux  (ITB49) 
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Figure  27.  Performance  Variation  of  the  Bundle  plus  Simulation 
Heaters  at  a  10  cm  Pool  Height  During  a  Decreasing  Heat  Flux  (DTB49) 
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Figure  28.  10  cm.  Performance  Variation  of  Tube  1  in  a  Pool  of  R- 
113  When  Influenced  by  a  Increasing  Number  of  Heated  Tubes  in  a 
Turbo-B  Bundle  During  a  Increasing  Heat  Flux 


63 


Heat  Flux/(W/m'‘2) 


100000 


10000 


1000 


100 


— 

B 

B 

■■1 

■II 

■taaaHHH 

IIMM 

— 

h. 

■ 

mtt 

■ 

1 

ai 

i 

iBMBani 

immm 

■ 

■ 

aaa 

III 

■1 

■ 

H 

■1: 

■ 

1 

1 

IMHHI 

B 

m 

■ 

III 

Uj 

■ 

■ 

Hi 

m 

m 

■ 

1 

9 

1 

III 

■■■1 

n 

■ 

1 

II 

iW 

H 

1 

1 

91 

1 

I 

1 

III 

_ AO/  r\ii 

■ 

1 

II 

z 

■ 

1 

1 

1 

1 

1 

III 

Ti 

jrbo- 

B  E 

un 

dk 

r 

IHI 

c 

3 

1 

1 

1 

mi 

mi 

1 

1 

1 

1 

II 

■ 

■ 

1 

1 

1 

1 

IM 

i 

1 

1 

1 

1 

1 

1 

1 

1 

II 

1 

1 

1 

1 

1 

!■ 

1 

1 

1 

1 

■ 

■ 

1 

1 

r 

1 

II 

■  1T(DTB04) 

•  2T(DTB05) 

■  3T(DTB06) 

•  4T(DTB07) 

■  5T(DTB27) 

D  B(DTB48) 

a  &i-S(DTB49) 

■■■■ 

■ 

■i 

ai 

a 

ai 

laaaaaa 

■■■■ 

II: 

■ 

1 

II 

!■■■■ 

■■■■■ 

_ 

H 

l|: 

■ 

1 

II 

!■■■■ 

■■■I 

■1 

■ 

■ 

l|a 

■ 

1 

1 

■ 

■ 

III 

■ 

1 

■ 

1 

II 

■ 

■ 

1 

■1 

■ 

1 

c 

ill 

■ipm 

iimH 

HI 

■ 

1 

1 

1 

■ 

1 

1 

1 

!■ 

■ 

1 

1 

1 

1 

■ 

1 

1 

11 

1 

1 

II 

mM 

■■ 

■■■■■■■■ 

■i 

91 

a 

Bl 

■■■■■a 

m 

■1 

If  ■■■M 

■ 

1 

ti 

n 

m 

■I 

liMMi 

■ 

1 

1 

iMBai 

H 

■ 

11 

mil 

■ 

1 

1 

1 

■ 

1 

1 

■ 

1 

1 

1 

■ 

■ 

1 

!■■■ 

■ 

■ 

1 

1 

■ 

1 

1 

IHI 

B 

■ 

1 

1 

Bl 

n 

■ 

1 

1 

1 

im 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

II 

1 

1 

1 

1 

1 

II 

1 

H 

1 

i 

.1 


10 


100 


(Twall-Tsat)/(K) 


Figure  29.  10  cm.  Performance  Variation  of  Tube  1  in  a  Pool  of  R> 
113  When  Influence  by  a  Increasing  Number  of  Heated  Tubes  in  a 
TurbO'B  Bundle  During  a  Decreasing  Heat  Flux 
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Figure  30.  Performance  of  Tube  1  at  a  20  cm  Pool  Height 
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Figure  33.  Performance  Variation  of  Tubes  1  to  3  at  a  20  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB35) 
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Figure  34.  Performance  Variation  of  Tubes  1  to  3  at  a  20  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB35) 
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Figure  37.  Performance  Variation  of  Tubes  1  to  5  at  a  20  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB37) 
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Figure  38.  Performance  Variation  of  Tubes  1  to  5  at  a  20  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB37) 
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Figure  40.  Performance  Variation  of  the  Bundle  at  a  20  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB38) 
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Figure  41.  Performance  Variation  of  the  Bundle  plus  Simulation 
Heaters  at  a  20  cm  Pool  Height  During  a  Increasing  Heat  Flux  (1TB40) 
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Figure  46.  Performance  Variation  of  Tubes  1&2  at  a  0  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB52) 
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Figure  48.  Performance  Variation  of  Tubes  1  to  3  at  a  0  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB53) 
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Figure  49.  Performance  Variation  of  Tubes  1  to  3  at  a  0  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB53) 
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Figure  50.  Performance  Variation  of  Tubes  1  to  4  at  a  0  cm  Pool 
Height  During  a  Increasing  Heat  Flux  (ITB55) 
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Figure  53.  Performance  Variation  of  Tubes  1  to  5  at  a  0  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB56) 
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Figure  55.  Performance  Variation  of  the  Bundle  at  a  0  cm  Pool 
Height  During  a  Decreasing  Heat  Flux  (DTB57) 
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Figure  56.  Performance  Variation  of  the  Bundle  plus  Simulation 
Heaters  at  a  0  cm  Pool  Height  During  a  Increasing  Heat  Flux  (ITB58) 
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Figure  57.  Performance  Variation  of  the  Bundle  plus  Simulation 
Heaters  at  a  0  cm  Pool  Height  During  a  Decreasing  Heat  Flux  (DTB58) 


92 


Heat  Flux/{W/m''2) 


Figure  58.  0  cm.  Performance  Variation  of  Tube  1  in  a  Pool  of  R- 
113  When  Influenced  by  an  Increasing  Number  of  Heated  Tubes  in  a 
Turbo-B  Bundle  During  a  Increasing  Heat  Flux 
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Figure  62.  Performance  Variation  of  Tubes  1  to  3  at  a  10  cm 
Pool  Height  During  a  Decreasing  Heat  Flux  with  Auxiliary 
Heaters  off  (DTB06T) 
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Figure  65.  Nucleation  Site  Deactivation  with  a  1  kW  Auxiliary 
Heater  Setting  and  a  20  cm  Pool  Height 
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Figure  67.  Nucleation  Site  Deactivation  with  a  1  kW  Auxiliary 
Heater  Setting  and  a  10  cm  Pool  Height 
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Figure  68.  Nucleatlon  Site  Deactivation  with  a  1  kW  Auxiliary 
Heater  Setting  and  a  0  cm  Pool  Height 
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Figure  70.  Nucleation  Site  Deactivation  with  a  1  kW  Auxiliary 
Heater  Setting  and  a  0  cm  Pool  Height 
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Figure  71.  Nucleation  Site  Deactivation  with  a  3  kW  Auxiliary 
Heater  Setting  and  a  10  cm  Pool  Height 
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Figure  73.  Nucleation  Site  Deactivation  with  a  3  kW  Auxiliary 
Heater  setting  and  a  10  cm  Pool  Height 


Heat  Flux/(W/m*2) 


100000 


10000 


1000 


100 


“n 

“1 

■ 

■ 

■■ 

M 

■ 

i4l 

IMHHH 

F— 

■■ 

■1 

■ 

■ 

■II 

H 

HI 

■ 

■ 

ia 

IMHHHi 

_ 1 

Ml 

M 

B 

■I 

IHHMHH 

M 

M 

■ai 

■ 

■ 

■ 

II 

immm 

“j 

■ 

■ 

B 

II 

■ 

1 

i 

■ 

1 

II 

!■■■ 

H 

m 

■ 

1 

1 

1 

■ 

8 

if 

lllll 

III 

m 

1 

1 

II 

■ 

U 

1 

1 

1 

IMHH 

m 

■ 

1 

1 

11 

R 

-113 

■ 

1 

1 

II 

■ 

1 

1 

1 

1 

IHHi 

H 

1 

1 

a 

li 

0%  Oll^ 

Turbo-B  E 

lun 

n 

dU 

II 

' 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

II 

1 

1 

||[| 

n 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

■ 

1 

1 

1 

1 

1 

1 

n 

1 

1 

1 

1 

1 

[[ 

1 

1 

1 

1 

1 

WM 

■ 

■ 

■1 

L_ 

M 

■ 

■1 

■■■HMi 

■i 

M 

M 

■II 

■ 

■ 

■ 

II 

■ 

■ 

■a 

IHHHHi 

m 

■ 

m 

iia 

■ 

■ 

II 

■ 

■1 

IMHHi 

M 

ai 

III 

— 

■ 

1 

II 

■ 

1 

1 

■ 

1 

III 

m 

■ 

1 

1 

1 

IHMi 

Ml 

m 

■ 

1 

1 

1 

■ 

1 

1 

III 

m 

■ 

1 

1 

1 

IH 

m 

1 

1 

1 

1 

a 

a 

III 

■ 

1 

1 

1 

1 

■■ 

H 

M 

1 

1 

1 

] 

1 

1 

1 

III 

■ 

1 

1 

1 

■ 

1 

1 

1 

n 

o  ios(rrB22) 

■  20s(ITB23) 

A  1  min(rTB24) 

1 

1 

1 

i 

1 

1 

1 

1 

1 

^IHHi 

■■■■■■ 

■ 

•h 

LJMana 

H 

■■ 

aasHM 

HI 

■ 

1 

II 

!■■■■ 

__j 

■I 

Ml 

M 

■ 

■1 

laHHi 

■ 

■ 

II 

■ 

1 

■1 

Mmm 

r 

n 

■ 

1 

II 

■ 

1 

■ 

1 

■ 

■ 

1 

1 

H 

1 

1 

1 

1 

■ 

1 

■ 

■ 

■ 

1 

■ 

1 

1 

■ 

■ 

a 

a 

■ 

■ 

— 

■ 

■ 

■ 

1 

1 

1 

WM 

■ 

1 

1 

llll 

1 

1 

1 

1 

lllll 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

.1 


10 


100 


(Twall-Tsat)/(K) 

Figure  74.  Nucleation  Site  Deactivation  with  a  3  kW  Auxiliary 
Heater  Setting  and  a  10  cm  Pool  Height 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

1.  When  secured  overnight,  the  bundle  exhibited  significant  temperature  overshoot 
for  increasing  heat  flux  (approx.  10*’C).  The  hysteresis  loop  formed  by  increasing 
followed  by  decreasing  heat  flux  was  very  repeatable. 

2.  In  the  natural  convection  region,  activating  lower  tubes  within  the  bundle  did  not 
affect  the  upper  tubes.  However,  the  heat-tran^er  coefficient  decreased  as  one  moved 
down  the  bundle. 

3.  Activation  of  the  lower  tubes  caused  the  upper  tubes  to  nucleate  earlier  (i.e.  at  a 
lower  heat  flux).  As  each  tube  successively  nucleated,  the  heat-transfer  performance  of  the 
remaining  lower  tubes  irtqnoved. 

4.  In  the  nucleate  boiling  region,  the  poformance  of  the  top  tube  was  the  best  This 
perframance  dropped  off  steadily  as  one  moved  down  the  bundle. 

5.  Variation  in  liquid  height  above  the  bundle  (finom  0  to  20  cm)  significantly  delayed 
die  point  of  incipience. 

6.  For  zero  pool  subcooling,  nucleation  site  deactivation  was  greatly  dependent  on 
both  pool  height  (i.e.  pressure  head)  and  tube  turn  off  time.  At  zero  pool  height, 
deactivation  was  negligible  for  all  time  whoeas  for  10  and  20  cm,  there  was  complete 
deactivation  in  60  and  15  seconds  respectively  (i.e.  decreasing  with  increasing  pool 
height). 

7.  At  a  typical  heat  flux  of  20  kW/m^,  the  Turbo-B  tube  bundle  heat-transfer 
coefficient  was  qrproximately  four  times  that  of  a  smooth  tube  bundle. 
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8.  The  uncertainty  of  the  data  was  greatest  at  low  heat  flux  and  low  wall  superheat 
due  to  the  measurement  being  of  the  order  of  accuracy  of  the  themMxouples. 

B .  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  Due  to  the  way  a  real  flooded  evaporator  operates,  it  would  be  preferable  to  heat 
the  tubes  within  the  bimdle  with  hot  water. 

2.  The  pool  height  tests  should  be  repeated,  but  by  varying  the  vapor  pressure  above 
the  liquid,  the  local  pressure  at  the  tube  surface  should  be  maintained  at  a  constant  value. 

3.  A  pressure  test  should  be  conducted  and  a  pressure  relief  valve  fitted  to  enable  the 
testing  of  alternative  refrigerants  such  as  R-124  which  has  a  higher  vapor  pressure. 

4.  Attention  needs  to  be  given  to  the  question  of  refrigerant  disposal.  There  are 
reclamation  projects  undertaken  by  most  manufacturers;  however,  a  method  still  needs  to 
be  found  to  remove  the  refrigerant  from  the  apparatus  into  a  container  for  reclamation. 

5.  Neutrally  buoyant  particles  should  be  placed  in  the  pool  to  facilitate  study  of  the 
circulatitMi  patterns  witiiin  the  bundle  in  more  detail. 

6.  A  high  speed  camera  should  be  used  to  study  the  nucleation  process  and 
circulation  patterns  in  more  detail  (the  mechanical  engineering  department  is  hoping  to 
purchase  such  a  camera  in  the  near  future). 

7.  A  sight  glass  should  be  fitted  to  the  refrigerant  storage  tank. 
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APPENDIX  A:  LISTING  OF  DATA  FILES 
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Number  of  Pairs  of  Number  of  Auxiliary  Pool  Height  Time  Tube 
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IkW 


Number  of  Pairs  of  Number  of  Auxifiary  Pool  Height  Time  Tube 
File  Name  Data  Points  Instrumented  Active  Tube  Simulation  Heater  Setting  in  centimeters  Heaters  wer 


(1.611 -0.0031  *(T))*1 00 


APPENDIX  B;  PROPERTIES  OF  R-113 
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APPENDIX  C:  SAMPLE  CALCULATIONS 


Data  set  number  1  of  run  DTB03  (Decreasing  heat  flux,  Turbo-B  tube.  Surface 
preparation  D)  was  used  for  the  sample  calculations  in  order  to  validate  the  program  used 
for  data  acquisition  DRP4A.  The  working  fluid  was  R-1 13. 

1.  Test  tube  dimensions 

=  11.60  mm 

Djj  *  14.15  mm 
Dj  =  12.70  mm 
L  =203.2  mm 
Lu  =  25.4  mm 

2.  Measured  Parameters 
Tl=  54.41  ®C 

T2  =  55.25  ®C 
T3  =  55.80  ®C 
T4  =  55.15  ®C 
T5  =  55.04®C 
T6  =  55.38®C 
Tldl=47.46*C 
Tld2  =  47.48®C 
Aas  =  3.838V 
Vas  =  3.958V 
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3.  Calculations 


The  heater  power  is  first  calculated  firom 

q=Vas(V)  X  Aas(V)  x  60(V/V)  x  1(A/V) 

Note:  The  multiplication  factors  of  volt  and  amp  sensors  are  60  and  1,  respectively. 


Therefore: 

q=  3.958V(3.838V)(60V/V)(1A/V) 
=  911.45W 


The  tube  inside  wall  temperature  is  obtained  from  the  average  of  all  six  thermocouple 
readings. 

_  6 
Twi  =iy  Tn 
^n-l 

*  1/6  (54.41  +  5525  +  55.80  +  55.15  +  55.04  +  55.38) 

=  55.17‘’C 


The  tube  outside  teiiq)erature  is  calculated  by  knowing  tiie  inside  wall  temperature  and 
using  Fouriers  Conduction  law .  Unifonn  radial  conductitm  is  assumed. 


Two  =  Twi  - 


qD«(^)l 

_ ^tc 

2x(kcu)(L) 
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where  the  second  term  on  the  right  hand  side  is  the  Fourier  conduction  term.  If  we 
define  this  term  as 


and 


0  = 


27l(kcu)(L) 


Sb  =  Two  -  Tsati, 


where  kcu  is  the  thermal  conductivity  and  is  calculated  as  follow 
kcu  =  434.0  -  [0.1(Twi-K)l 
=  434.0 -[0.1(328.32)1 
=  401.17  WAnK 


Now, 


_  _  911.45an(|f^)l 

Two  =  Twi - 


2ii(401.17)(0.2032) 
=  (55.17 -0.3567)®C 
=  54.82«C 


The  liquid  saturation  teniq>erature  at  the  top  of  the  tube  bundle  is 

Tsat='n^^  tTId2 

2 

2 

=  47.47®C 
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In  order  to  calculate  the  local  saturation  ten^>erature  for  each  tube,  coirection  factors 
are  needed  to  account  for  hydrostatic  pressure  differences  between  the  tube  locations  and 
the  liquid  free  surface.  This  difference  is  calculated  by: 

AP  =  p(g)(ht) 

where  the  density  is  calculated  by 

p  =  1620.7479  -  Tsat(2.2186346  +  Tsat(2.3578291E-3)) 

Therefore,  for  Tsat  =  47.47®C 
p=  1510.1  kg/in3 

So,  for  the  top  tube  in  the  bundle  which  is  0.0124  m  below  the  thermocouples 
measuring  the  pool  temperature, 

dP=  1510.1  (9.81)(0.0124) 

AP=  183.69  Pa 

For  a  183.69  Pa  pressure  difference,the  corrected  saturation  temperature  is  obtained 
by  adding  0.054**C  (obtained  from  standard  tables  for  R-113  or  by  equations  shown  in 
Appendix  F.)  to  Tsat  The  corrected  Tsat  is: 

Tsatc  *  (Tsat  +  0.054)®C 
=  (47.47 +  0.054)«C 
Tsatc  *47.52*0 
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Therefore,  the  wall  superheat  can  be  obtained  by  the  following 


0b  =  Two  -  Tsatc 
=  (54.82  -  47.52)®C 


=  7.30®C 


Now  that  the  wall  superheat  is  known,  we  need  to  calculate  the  heat  flux  and  the  heat 
transfer  coefficient  To  do  this,  we  know  that  die  tube  is  12  inches  long  and  is  heated  in  a 
ei^t  inch  center  porticm  of  the  tube.  The  unheated  lengths  of  the  tube  are  a  one  inch  and  a 
three  inch  section  on  oppoate  ends  of  the  tube.  Hiese  unheated  lengths  have  a  fin  effect 
during  the  heat  transfer  process  to  the  evaporating  refrigerant  In  order  to  account  for  this, 
the  following  procedure  was  adopted  for  both  the  one  and  three  inch  sections.  Calculatims 
are  shown  below  for  the  one  inch  section . 

Heat  transfer  from  die  unheated  end  is  calculated  as  heat  from  the  base  of  a  fin: 

qf*(hbXpxkcuXAc)®’^x  6bXtanh(nxLc) 

where 

p  =  icxDo 
=  icx  14.1SE-3m 
=  44.45E-3m 

now 

Ac  =  ?  (0.0141?  -  0.01270>) 

4 

Ac  =  3.0578E-5m2 
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The  conected  length  of  unenhanced  surface  at  the  end  was  calculated  as  follows 

1,  =  0.0254 

Lc  =  0.02S8  m 


hb  is  the  natural  convection  heat  transfer  coefficient  of  the  fin  like  ends  and  was 
calculated  by  using  the  Churchill-Chu  [Ref.  24]  correlation  for  natural  convection  from  a 
smooth  horizontal  cylinder,  as  modified  by  Pulido  [Ref.  27]. 


gxPxE^xObX  tanh(n  xLc) 

1 

6 

h.  -  Jl 

n  A  4.0  V 

vxaxLcxn 

\/«0 

[i  .[0^]^- 

- 

2 


where 


n  = 


hbXp 


Icqi  xA^ 


0.5 


Therefore  an  iterative  technique  was  necessary  to  calculate  hb.  The  iterative  technique 

JK_ 

used  was  to  assume  hb  was  190  j„2i(  and  to  continue  the  iteration  until  successive  values 
are  within  0.001  of  each  other.  The  fluid  physical  properties  are  calculated  at  the  vapor 
mean  film  temperature,  given  by  the  following  equation. 
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Tfita  =  51.17®C=324.32®K 

Fot  R-113,  the  physical  properties  are  given  in  the  program  by 

Dynamic  viscosity,  Tfiim  ® 

H=1.34E-5xlo[T«”£ii] 

\i  =  4.88E-4  kg/m-s 
Specific  heat,  Tnim  in  ®C 

Cp  =  929.0  +  (1.03  X  Tfite) 

Cp  =  981.71  J/kg-K 

Density,  Tfiim  in  ®C 

p  =  1620.7479  -  Tfii„f2.2186346  +  Tfiim(0.0023578291)] 
p  1501.1  kgAiP 

Thermal  conductivity  of  R-113,  Tfjm  in  ®C 

k  *  0.071  -  (0.000261  x  Tgim) 

k  =  5.765E-2W/mK 

Prandtl  number 

Pr  =  q,xt 

Pr  SB  981  71  X  4i88E-4- 
Pr  =  8.31 
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Thermal  mpansion 

p  =  -(l/p)(Ap/AT) 


t)5i.07  =  1501.29  kg/m^ 
p  =  -(l/1501.1)(0.49/0.2) 
p=  1.63E-3  (1/K) 
Kinonadc  viscosity 


P 

1501.1 

v  =  3.25E-7m2/s 


Thermal  difiusivity 

a= — 
pxCp 


5.7645E-2 
”  1501.1  X  981.71 


a  =  3.912E-8m2/s 


PsiJt?  =  1500.8  kg/m^ 


Knowing  the  above  properties,  the  heat  transfer  coefficient,  hb,  can  be  obtained  by 
iteration. 

hb  =  185.42  W/m2« 


knowing  this  we  calculate  n 


hbXP 
_keu  Ac 


0.5 
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[  185.42  X  44.45E-31 
"[401.17  X  3.0578E-5 

n  =  25.92 

then  we  can  obtain  qf 

qj s  (hb  xp  xkcu X Ac)®*^  x  ©b  x tanh(n xLc) 

qf  =  (185.42  X  44.45E-3  x  401.17  x  3.0578E-5)®-^  x  7.3  x  tanh  (25.92  x  0.0258) 
qf=1.36W 

The  corresponding  results  for  the  three  inch  section  are 
hb  =  158.27  W/m2.K 
qf  =  2.04W 

Therefcxe  the  heat  tnuisfer  through  the  heated  lengA  of  the  tube  is 
qs  -  q  -  qf(l  inch  section)  -  qf(3  inch  section) 
qs  *(91 1.45 -1.36  -  2.04)  W 
qs  =  908W 

and  the  heat  flux  and  the  heat  transfer  coefficient  are  as  follows 

q"  *  qs/As 

=  qi/(JCxDoXL) 

*  908/(JC  X  0.01415  X  0,2032) 

-908/9.033E-3 

*1.005E+5W/m2 
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and  finally  the  heat  transfer  coefficient 


A,  X  (Two-Tsat) 

=  908/(9.033E-3x7.3) 
h  =  1.376E+4  W/m2.K 


APPENDIX  D:  UNCERTAINTY  ANALYSIS 


The  same  data  run  (DTB03)  was  chosen  for  the  uncertainty  analysis.  Therefore,  the 
measured  and  calculated  parameters  found  in  the  sample  calculation  were  used  in  this 
section.  The  uncertainty  analysis  performed  was  for  a  high  heat  flux  but  the  same 
procedure  could  be  performed  at  any  heat  flux  to  determine  the  uncertainty  bands.  All 
uncertainties  are  presented  as  a  percentage  of  the  calculated  parameter.  The  uncertainty 
associated  with  the  experimental  parameters  is  calculated  from  the  equation  suggested  by 
Kline  and  McQintock  [Ref  28].  For  example,  if 

R  *  R  (xj,  X2,  •  • . , 

dien 

SR=  (^.)  + 

oXi  0X2 

where: 

$R  s  uncertainty  of  the  desired  dependant  variables 

Xn  =  measured  variables 

5xn  =  uncertainty  in  measured  variables 

The  boiling  heat  transfer  coefficient  is  given  by 

h - - 

A,  x(Two-Tsat) 

where 

_  _  q[i”(fe)i 

Two-Twi-;— -W- 

27c(kcu)(L) 
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In  the  above  equation,  the  second  term  on  the  right  hand  side  is  usually  called  the 
Fourier  heat  transfer  conduction  term  If  we  deHne  this  as  O,  then 

27iQtcu)(L) 

and 

0b  =  Two  -  Tsatc 

With  this  notation,  the  uncertainty  in  the  heat  transfer  coefficient  is  obtained  using  the 
following  equation. 

5k = %  [5^1  \  f  ^1  %  f  Ml  %  f  SImI  2 

h  [LqJ  La,J  [  J  [e^J  [  Gb  J  . 

where 

q  =  VxI 

q  *  V(V)  X  I(V)  X  60(V/V)  X  1(A/V) 

and  the  uncertainty  is 

i-ifi’-ffli- 

The  accuracy  in  dte  voltage  and  current  sensors  are  as  follows 

5Vas  =  ±0.05  V  5Aas  =  ±0.025  A 
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FrcHn  the  sample  calculation  section 
Vas  =  3.958V 


Aas  =  3.838V 


Therefore, 

8q/q  « ((8Vas/Vas)^  +  (8Aas/Aas)^) 
8q/q  =  ((0.05/3.958)2  +  (0.025/3.838)2) 


8q/q  =  1.42  percent 


Calculation  of  the  surface  area  and  die  uncertainty  of  it  are  done  as  follows 
As  =  icxDoXL 


Knowing  the  dimensions  of  the  tube  firom  the  manufacturer  and  estimated 
inaccuracies  from  work  shop  tools  and  human  error,  the  uncertainty  was  calculated. 

Dimensions 

Do  « 14.15  mm  L  »  203.2  mm 

Inaccuracies  in  measurements 

SDo«0.1mm  8L*0.2mm 
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Uncertainty  analysis  performed 

5A,/A,  =  ((5Do/Do)*  +  (JSULf)^^ 
bAJA,  =  ((0.1/14.15)2  +  (0.2/203.2)2)®-^ 
bAJAt  =  0.7135  percent 

The  uncertainty  calculation  for  the  Fourier  conduction  term  is  given  below 


kcu  was  calculated  using 

kcu  =  434.0 -[0.1(Twi-K)] 

=  434.0 -[0.1(328.32)1 
=  401.17  W/mK 
and  its  uncertainty 

5kcu=[(0.1  x5Twi-K)*]2 

^wi  and  CTsat  ate  obtained  using  uncertainties  in  the  thermocouple  readings. 
Average  wall  inside  ten^erature  Twi  was  obtained  taking  the  average  of  six  thermocouple 
readings  inside  the  tube.  The  uncertainty  associated  with  this  variable  is 


where  STc  was  obtained  as  ±  0.5*’C  from  "Omega  1987  Complete  Temperature 
Measurements  Handbook  pp.  T-37.  Using  engineering  judgement  this  seems  high  will  use 
0.1®C. 
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8Twi= 


5Twi  =  0.04«C 

Saturation  temperature  was  obtained  by  taking  the  average  of  two  thermocouple 
readings  and  the  uncertainty  in  this  temperature  was  calculated  fix>m  the  following  equation. 


STsat  =  0.07®C 


Knowing  the  uncertainty  in  the  temperatures,  we  can  now  calculate  the  uncmainties 
in  the  following. 

6kcu=[(0.1  x^wi-K)*]2 
5kcu=[(0.1  X  273.19)^]  2 
5kcu  =  27.32  W/inK 


Now  we  can  calculate  the  uncertainty  in  the  Fourier  conductitxi  tenn 
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^  =  6.96  percent 

0 

frun  the  sample  calculations  we  know  that  O  =  0.3567  or 
50  =  0.3567  X  6.96  percent = 0.025  ®C 


thoefore: 


[8a 

2 

+ 

8A. 

2 

+ 

[STai] 

2 

+ 

50 

2 

+ 

[£ISaL] 

2 

LqJ 

.  . 

Ob 

[0bJ 

0b 

^  -  [[0.0142]  >  +  [0.07 135]  2 + [£^]  * + [0^]  *  +  [£^]  *]  2 

^  =  7.4  percent 
h 


Table  3  shows  the  results  of  the  uncertainty  analysis  perfoimed.  The  high  and  low 
heat  flux  correspond  to  the  approximate  values  of  1.005E+5  W/m^  and  597  W/m^ 
respectively. 
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TABLE  3.  UNCERTAINTY  ANALYSIS  RESULTS 


Variable 

High  Heat  Flux 

Low  Heat  Flux 

AT 

7.3 

0.38 

Twi 

55.17 

48.04 

Tsat 

47.47 

47.6 

SY 

V 

1.26% 

16.4% 

SL 

I 

0.65% 

8.3% 

8q 

q 

1.42% 

18.% 

8Do 

Do 

0.707% 

0.707% 

SL 

L 

0.098% 

0.098% 

5A, 

A. 

0.714% 

0.714% 

kcu 

6.8% 

6.8% 

o 

6.96% 

21.9% 

StL 

h 

7.4% 

22.4% 
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APPENDIX  E:  OPERATING  PROCEDURES 


SYSTEM  LIGHT-OFF 

1 .  Power  to  the  8  ton  refrigeration  unit  is  provided  by  the  breakers  located  in 
the  main  distribution  panel  located  in  the  laboratory.  These  breakers  were  never  secured. 
However,  if  power  to  this  panel  was  lost  then  these  breakers  must  be  reset 

2.  Turn  the  switch  on  the  refirigeration  unit  control  panel,  located  in  front  of  the 
refrigeration  unit  to  the  "auto"  position  after  passing  through  "on"  position.  This  switch  is 
always  left  on,  unless  unit  was  taken  down  fOT  long  repairs. 

2.  Push  the  start  button  in  the  control  box  for  the  recirculation  punq>.  This 
control  box  is  located  on  the  bulkhead  above  the  recirculation  pump  in  the  outside  area 
adjacent  to  the  refrigeratitxi  unit 

3.  Set  the  desired  temperature  on  the  roughly  graduated  Fahrenheit  scale  on  the 
control  panel  thermostat  It  requires  approximately  one  hour  to  chiU  the  sump  S°C  The 
thermometer  of  the  ethylene  glycol  sump  must  be  monitored  to  ensure  the  desired  sump 
temperature  is  attained  and  maintained.  Slight  adjustments  in  the  refrigeration  unit 
tiiermostat  can  be  expected  due  to  the  coarseness  of  its  scale 

4.  Energize  the  desired  pumps  by  switching  on  the  breakers  in  the  main 
distribution  power  panel  Once  the  power  is  energized  in  the  main  power  panel  turn  on  the 
punq)  motors  by  pushing  down  on  the  arm  of  the  appropriate  breaker  box  for  the  pumps 
located  on  the  bulkhead  next  to  the  ethylene  glycol  sump.  The  pumps  are  marked 
"auxiliary  condenser"(Pump  #2)  and  "instrumented  tube  condenser"(Pump  #1), 
respectively  in  the  breaker  box.  Flow  in  the  auxiliary  cmdensate  system  can  be  controlled 
with  the  individual  globe  valves  located  at  the  coil  penetrations  on  the  apparatus.  The 
auxiliary  condenser  will  produce  the  fastest  adjustments  to  system  pressure  if  pressure  is 
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rising  too  quickly.  Row  through  the  instrumented  tubes  can  be  controlled  by  the  ball 
valves  located  at  the  bottom  of  the  respective  flow  meter. 


5 .  Energize  the  heater  variacs  desired  by  switching  on  the  breakers  (Auxiliary, 
Simulation,  and/or  Bundle)  in  the  main  distribution  power  panel  (located  near  entrance  of 
lab)  and  individual  breakers  for  each  of  these  in  the  power  distribution  box  (near 
apparatus). 

6.  After  ensuring  that  the  breakers  fcv  the  heater  tubes  desired  are  in  the  "on” 
position,  adjust  variac/s  to  the  desired  amount  of  heat  flux.  Monitor  apparatus  pressure 
through  the  pressure  gage  at  the  top  of  the  apparatus,  ensuring  system  pressure  does  not 
exceed  30  psig. 

SYSTEM  SHUTDOWN 

1 .  Turn  all  variacs  to  the  zero  position  and  switch  off  all  breakers  in  the  power 

panels. 

2.  If  apparatus  will  not  be  operas  for  an  extended  period,  turn  the  switch  on 
the  refrigeration  control  panel  to  the  "off'  position  afro’ passing  through  "on”. 

3.  Allow  the  recirculation  pump  to  operate  for  at  least  five  minutes  after 
switching  off  the  refrigeration  unit  to  dissipate  any  back  pressure  in  the  system. 

4.  Turn  the  breakers  for  the  pumps  to  the  off  position  at  the  switch  boxes  near 
the  ethylene  glycol  sump,  and  then  secure  powo*  at  the  main  distribution  power  panel. 

EMERGEMCY-  SHUTPaWN 

1.  Secure  all  power  at  the  main  distribution  power  panel 

2.  Pull  fire  alarm 

3.  Evacuate  building 
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APPENDIX  F:  PROGRAM  DRP4A 

1000!  FILE  NAME=  0RP4A 

10041  DATE:  November  22,  1988 

10081  REVISED:  Nay  1991  (S.  MEMORY) 

10121 

1015  BEEP 
1020  PRINTER  IS  1 
1024  ldp*0 
1028* 

1032  PRINT  USING  “4X, ""Select  option  default  is  0:""" 

1035  PRINT  USING  "SX,“"0  Taking  data  or  re-processing  previous  data 

1040  PRINT  USING  "GX,“"1  Plotting  data  on  Log-Log 

1044  PRINT  USING  "SX,"“2  Plotting  data  on  Linear""" 

1048  PRINT  USING  “GX,""3  Purge""" 

1052  PRINT  USING  ■■GX,‘"4  Fixup""" 

1055  PRINT  USING  "GX,""5  Move""" 

1050  PRINT  USING  "GX.'-G  Comb""" 

1054  PRINT  USING  "SX,""7  Read  Plot""" 

10581 

1072!  IDP  IS  A  PROGRAM  VARIABLE  TO  SELECT  A  SUBROUTINE 
1075  INPUT  I  dp 

1080  IF  ldp>0  THEN  CALL  Main 
1084  IF  Idp-1  THEN  CALL  Plot 
1088  IF  I  dp-2  THEN  CALL  PI in 
1092  IF  Idp-3  THEN  CALL  Purg 
1095  IF  Idp-4  THEN  CALL  Fixup 
1100  IF  Idp-S  THEN  CALL  Move 
1104  IF  I dp-6  THEN  CALL  Comb 
1108  IF  Idp-7  THEN  CALL  Readplot 
1112  END 

1116! 

1120  SUB  Main 

1124!  ICAL-THERMOCOUPLE  CALIBRATION 
1128  COM  /Cc/  C(7) 

1132  DIM  Enf<35>,T(3S),01a<6 ),02a<6  >.0ia< 6 ),0oa< 6  >,La< 6  > .Lua( 6  > ,Kcua< 6  >,Et ( 19 > , 
Ldtc( 20 ), Volt ( 2  >, Anp< 11  ),Tua<5>,Tu(5  ),Theta(5 ),THetab( 5  > ,0<5  )  ,01( 5 ) ,0dp(  5 ) 

1136  DIM  Htube(5) 

1140! 

1144'  THERMOCOUPLE  ARRAY  <C(  ))  INITIALIZATION 
1148  DATA  0.10086091,25727.94369,-767345.8295,79025595.81 
1152  DATA  -9247486589, 6. 97688E+11, -2. 66192E+13. 3. 94078E+14 
1155  READ  CO) 

1150! 

11641  PRINT  HEADER  AND  INITIALIZE  TIME  CLOCK 
1168  PRINTER  IS  701 
1172  BEEP 

1175  INPUT  "ENTER  MONTH,  DATE  AND  TIME  ( MM- DD : HH : MM  SS ) " , Da t e$ 

11801  OUTPUT  DIRECTED  TO  DATA  AOUISITION  SYSTEM  'HP  3497A ) 


1184 

OUTPUT  709; "TO": Date* 

1188 

OUTPUT  709; "TD" 

1192 

ENTER  709; Date* 

1196 

PRINT 

1200 

PRINT  "  Month,  date  and  tine  : 

”  ;DateS 

1204 

PRINT 

1208 

PRINT  USING  ■10X,""NOTE:  Program  name  : 

DRP4"" 

1212 

BEEP 

12201  ON  IS  THE  VARIABLE  FOR  DISC  NUMBER  FOR  RECORD  KEEPING  ONLY 

1224  INPUT  “ENTER  DISK  NUMBER", Dn 

1228  PRINT  USING  ■lGX,““Disk  number  -  ‘“.ZZ-iDn 

1232  BEEP 

1236  In-0 

1240  INPUT  "ENTER  INPUT  MODE  (0-3497A, 1-FlLE  )  0-DEFAULT" , In 
1244! 

1248!  INPUT  MODE  ZERO  IS  FROM  THE  DATA  AQUISITION  SYSTEM 
1252  IF  In-0  THEN 
1256  BEEP 

1260  INPUT  "6IUE  A  r AME  FOR  THE  RAW  DATA  FILE* ,D2f ileS 

1264  PRINT  USING  "lEX.^File  nane--  ,  14A“  !02f  i le* 

1268! 

1272!  CREATE  BOAT  FILE  ON  THE  MASS  STORAGE  MEDIA 
1276  CREATE  BDAT  D2 filet, 60 

1280!  CREATE  AN  INPUT/OUTPUT  LINK  TO  OPEN  FILES 
1284  ASSIGN  0File2  TO  02fileS 

1288! 

1292!  CREATE  DUMMY  FILE  UNTIL  Nrun  KNOWN 
1296  Dlfilet-'DUMMY" 

1300  CREATE  BOAT  01  filet, 60 

1304  ASSIGN  0Filel  TO  Olfilet 

1308  OUTPUT  OFileliDatet 

13121 

1316!  CREATE  A  PLOT  FILE 

1320  BEEP 

1324  INPUT  "GIUE  A  NAME  FOR  THE  PLOT  FILE’,Pfilet 

1328  CREATE  BDAT  Pfilet.30 

1332  ASSIGN  OPlot  TO  Pfilet 

1336  BEEP 

1340! 

1344!  IDTC  -  NUMBER  (TOTAL)  OF  DEFECTIVE  THERMOCOUPLES 
1348  INPUT  “ENTER  NUMBER  OF  DEFECTIVE  TCS  ( 0-OEFAULT >“ , Idtc 
1352!  LDTC  -  LOCATION  OF  DEFECTIVE  THERMOCOUPLE 

1356! 

1360  IF  Idtc-0  THEN 

1364  PRINT  USING  •16X,“*No  defective  TCs  exist*” 

1368  ELSE 

1372  PRINT  USING  *16X, * “Defective  THernocouples  Indicated  by  -99.99'““ 

1376  END  IF 

1380! 

1384  BEEP 

1388!  DEFECTIVE  THERMOCOUPLES  MAY  BE  IN  CHANNELS  40-69 

1392!  THEMOCOUPLES  ARE  ENTERED  AS  DEFECTIVE  BY  COMPUTER  CHANNEL  NO. 

1396!  JOTC-COUNTER  IN  LOOP  FOR  DEFECTIVE  THERMOCOUPLES 

1400! 

1404  IF  ldtc>0  THEN 

1408  FOR  Jdtc-0  TO  Idtc-1 

1412  INPUT  “ENTER  DEFECTIVE  TC  LOCATION  (BY  COMPUTER  CHANNEL  NUMBER) 

",Ldtc< Jdtc  ) 

1416  BEEP 

1420  NEXT  Jdtc 

1424  END  IF 

1428  PRINTER  IS  701 

1432  OUTPUT  0Kilel:Ldtc(«  ) 

1436  '  ^ 

1440  !  In-1  option  (THIS  OPTION  ALLOWS  DATA  ENTRY  WITH  DATA  FILE) 

1444  ELSE 
1448  BEEP 

1452  INPUT  “GIVE  THE  NAME  OF  THE  EXISTING  DATA  FILE" .02f ilet 
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1456  PRINT  USIN6  '16X,’"File  name-  " “ . 14A“ ; D2f i le* 

1460  ASSIGN  9File2  TO  □2file$ 

1464  ENTER  eFile2iNrun 

1468  ENTER  0File2 i DoldS,Ldtc( •  ) , Itt ,Bop ,Nht ,Natp ,Nrt ,Corr 
1472  BEEP 

1476  INPUT  "GIUE  A  NAME  FOR  PLOT  FILE“.Pfile$ 

1480  CREATE  BOAT  Pfilet,30 

1484  ASSIGN  SPlot  TO  PfileS 

1488  PRINT  USING  ■16X,""This  data  set  taken  on  :  ■ “ , UA" ; Doldt 

1432  BEEP 

1496  PRINTER  IS  1 

1500  PRINT  USING  MX,  “'SELECT  TUBE  TYPE*” 

1504  PRINT  USING  “6X,"0  SMOOTH'" 

1508  PRINT  USING  “6X,““1  FINNEOl 19/IN >  ”* 

1512  PRINT  USING  “6X,"2  HIGH  FLUX 

1516  PRINT  USING  ■6X,”3  TURBO-B  ”* 

1520  INPUT  lit 

1524  END  IF 

1528  IF  In-l  THEN  GOTO  1768 

1532  PRINTER  IS  1 

1536! 

1540  IF  Iin-0  TUEN 

1544  PRIN,  USING  ■4X,  "Select  tube  type"" 

1548  PRINT  USING  “6X,"  0  Smooth  “*“ 

1552  PRINT  USING  “BX,'*  1  FINNED  19/IN  (DEFAULT)"” 

1556  PRINT  USING  “6X,"  2  HIGH  FLUX**“ 

1560  PRINT  USING  “6X,"  3  TURBO-B"* 

1564  PRINT  USING  *6X,**  4  GROWTH"* 

1568  PRINT  USING  *6X,**  5  GROWTH"* 

1572  PRINT  USING  *6X,“*  6  GROWTH*" 

15761  ITT-TUBE  TYPE 

1580  INPUT  Itt 

1584  OUTPUT  •Fileliitt 

1568  END  IF 

1592  PRINTER  IS  701 

1596!  Itt-2 

1600  PRINT  USING  ■16X,"Tube  Type:  **,DD*:Itt 

1604! 

1608  BEEP 
1612  Bop*0 

1616  INPUT  "ENTER  BULK  OIL  X  <DEFAULT-0X>  *,Bop 
1620  OUTPUT  0Filel;Bop 

1624  PRINT  USING  *16X,"Bulk  oilX-*  '  ,DD*  iBop 
1628! 

1632  BEEP 

1636!  NHT-NUMBER  OF. HEATED  TUBES 
1640  Nht-S 

1644  INPUT  "Enter  number  of  heated  instrumented  tubes! defaul t“5 >" ,Nht 
1648  OUTPUT  SFneliNht 

1652  PRINT  USING  " IGX ,* "Number  of  heated  instrumented  tubes-" * ,00" :Nht 
1656  BEEP 
16601 

1664!  Natp-Number  of  active  dummy  pairs 
1668  Natp-0 

1672  INPUT  "Enter  number  of  active  dummy  pairs  (DefauIt-0 )” ,Natp 
1676  OUTPUT  SFileliNatp 

1680  PRINT  USING  "IGX, "Number  of  active  dummy  pairs-" ’ ,00" ; Natp 
1684  BEEP 
1688! 

1692!  NRT-NUMBER  OF  ADDED  HEATED  TUBES  TO  ENHANCE  BUNDLE  EFFECT 


1696  Nrt«0 

1700  INPUT  “Enter  number  of  added  heated  tubes  from  simulation  heaters(Default» 
0)*,Nrt 

1704  OUTPUT  OFileliNrt 

1708  PRINT  USING  “ 16X ,* “Number  of  added  heated  tubesCfrom  simulation  heaters 

“.DD'iNrt 

1712  BEEP 

1716! 

1720!  CORR  IS  CORRECTION  FOR  INSTRUMENTED  TUBE  HEIGHT 
1724  Corr-0 

1728  INPUT  "WrtNT  TO  CORRECT  TSflT  FOR  TUBE  HEIGHT  (0-VES<OEFflULT ), 1-NO )“ ,Corr 
1732  IF  Corr-0  THEN  PRINT  USING  "1GX,*"TSAT  is  corrected  instrumented  heat 

ed  tube  height  ““ 

1736  IF  Corr-1  THEN  PRINT  USING  “16X,‘“TSAT  is  NOT  corrected  for  instrumen 

ted  heated  tube  height*** 

1740  OUTPUT  0Filel:Corr 
1744  BEEP 

1748!  ILOU-INPUT  MODE:  LIQUID.  UAPOR.OR  LIQUID  UAPOR  AUERA6E 
17B2  Ilqv-0 

1756  INPUT  “SELECT  ( 0-LIQ( default ), l-UAP,2-( LIQ+UAP )/2 Ilqv 
1760! 

1764!  OlA-Oianeter  at  thermocouple  positions  < meters! 

1768  DATA  .0122,0.0098,0.0106.0.0116.0.0,0 
1772  READ  OlaO) 

1776  Dl-Dla(Itt) 

1780! 

1784!  02-0iameter  to  base  of  fins  (outside  die  for  smooth X meters > 

1788  DATA  .0158.0.0125,0.0158,0.01415.0,0,0 
1792  READ  D2a<*) 

1796  02-02a<Itt> 

1800! 

1804!  Di-Inside  diameter  of  unenhanced  ends  (meters) 

1808  DATA  .0132,0.0109,0.0116.8.0127,0,0.0 
1812  READ  Dla(*> 

1816  Oi-Oia(Itt) 

1820! 

1824!  Do-Outside  diameter  of  unenhanced  ends  (meters) 

1828  DATA  .015875,0.0125,0.015875.0.01415,0,0,0 
1832  READ  OoaO) 

1836  Do-Ooa(Itt) 

1840! 

1844!  L-Length  of  enhanced  surface  (meters) 

1848  DATA  .2032, .2032, .2032, .2032, .2032, .2032, .2032 
1852  READ  La(*) 

1856  L«La(Itt) 

18601 

1864'  Lu-CORRECTED  Length  of  unenhanced  surface  at  the  ends  (METERS) 

1868ILU-LFIN  +  THICKNESS/2 

1872  DATA  .0261 , .0254, .0264,0.0258,0,0,0 

1876  READ  LuaO) 

1880  Lu-Lua(Itt) 

1884! 

1888!  LU-corrected  length  of  3  inch  finned  like  end 
1892  DIM  Lva(6) 

1896  DATA  .0769, .0762, .0773,0.0766,0,0,0 
1900  RfAO  LvaO) 

1904  Lv-Lva(Itt) 

1908!  Kcua-Thermal  Conductivity  of  tube 
1912!  DATA  401,0,0,0,0,0,0 
19161  READ  Kcua(*) 
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1920'  Kcu-Kcua< Itt ) 

1924  A-PI»(0o*2-Di'2)/4 

1928  P-PI*Oo 

1932  J-1 

1936  Sx-e 

1940  Sy-0 

1944  Sxa-0 

1948  Sxy-e 

1952  Repeat:  ! 

1956! 

1960  IP  In-0  THEN 

19641  Dt Id-desired  tenperature  of  liquid 
1968  Otld-47.5  iR-113 

19691  Otld-2.2  IR-114 

1972  Ido-2 

1976  ON  KEY  0,15  RECOVER  1952 

1980  PRINTER  IS  1 

1984  PRINT  USING  “4X,‘ ‘SELECT  OPTION  ■“ 

1988  PRINT  USING  ‘SX.-'O-TAKE  DATA*** 

1992  PRINT  USIN6  *6X,**1-SET  HEAT  FLUX*** 

1996  PRINT  USING  *GX,**2-SET  Tsat  (DEFAULT  SET  FDR  R-113)““* 

2000  PRINT  USING  *4X,’*N0TE:  KEY  0  -  ESCAPE* ** 

20041  Ido-desired  option 

2008  BEEP 

2012  INPUT  Ido 

20161 

2020  BEEP 

20241  Set  default  value  for  input 

2028  IF  Ido>2  THEN  Ido-2 

20321  Take  data  option 

2036  IF  Ido-O  THEN  2440 

20401 

2044!  LOOP  TO  SET  HEAT  FLUX  (FOR  TOP  INSTRUMENTED  TUBE) 

2048  IF  Ido-1  THEN 

2052  Dqdp-1 00000 


2056 

PRINT  USING  •4X,“Qdp 

QDPsin 

Nrt 

Qdpaux 

Otof* 

2060 

PRINT  USING  •4X,“’(W/n*2) 

(W/n*2> 

(W/n*2) 

(U)*-* 

2064 

Err-1 

2068!  Reset, read  channel  25-30, autonatic  scaling 

2072!  Channel  25-aux  anps, 26-sin  anps,27-in8t  volts, 28-sin  volts, 29-aux 

volts, 30-inst  anps 


2076 

OUTPUT  709i*AR  AF25  AL30  VR5 

2080 

F:R  1-10  TO  11 

2084 

OUTPUT  709; "AS  SA* 

2088 

ENTER  709:Anp(I) 

2092 

NEXT  I 

2096 

FOR  1-0  TO  2 

2100 

OUTPUT  709; “AS  SA* 

2104 

ENTER  709; Volt!  I) 

2108 

NEXT  I 

2112 

OUTPUT  709: "AS  SA" 

2116 

ENTER  709:Anp(0) 

2120! 

Calculate  actual  heat  flux 

2124  ^ 

0(  0  )-60* Vol  t  (  0  )*Anp ( 0  ) 

2128 

Odp(0>-O(0)/(PI»D2»L) 

2132 

Osln-60*20*Volt( 1 >«Anp(ll > 

2136 

Qdpsin-Osin/ ( PI»D2* .2032»3  ) 

2140 

Oaux-60*20*Volt(2)*Anp( 10) 
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2144 

Odpaux-Qaux/ (  PI*. 0160*. 1776*4  ) 

2148 

Otot-Q( 0 )*Nht+05in+Qaux 

2152 

Nrt-0dpsin/0dp<0 ) 

2156 

IF  ABS(Aqdp-DQdp»Err  THEN 

2160 

IP  Aqdp>Dqdp  THEN 

2164 

BEEP  4000,. 2 

2168 

ELSE 

2172 

BEEP  250,. 2 

2176 

END  IP 

2188 

IF  Nrt<.l  THEN  Nrt-0 

2184 

IP  Odpaux<100  THEN  Qdpaux-0 

2188 

IP  Odpsin<100  THEN  Qdpsin-B 

2192 

PRINT  USIN6  •4X,2(MZ.3DE,2X).2X.(HDD.DD),2X,: 

(0  ),Odpain,Nrt 

.Qdpaux.Otot 

2196 

WAIT  2 

2200 

GOTO  2076 

2204 

END  IF 

2208 

END 

IF 

2212! 

22161 

LOOP 

TO  SET  Taat 

2220 

IP  Ido-2  THEN 

2224 

IF  Ikdt-1  THEN  2240 

2228 

BEEP 

2232 

INPUT  ‘ENTER  DESIRED  Tsat  (DEFAULT-47. 5  C  -  R- 

2233! 

INPUT  ‘ENTER  OFSIREO  Tsat  (DEFAULT-2.2  C  -  R- 

2236 

Ikdt-1 

2240 

Oldl-0 

2244 

Old2-0 

2248 

Nn-1 

2252 

Nrs-Nn  HOD  15 

2256 

Nn-Nn+1 

2260 

IF  Nrs-1  THEN 

2264 

PRINT  USING  UX,’*  DTsat  Tldl 

Tv«i 

Tvab  Tlav 

2268 

END  IF 

2272! 

Read  therMocoupla  voltages  for  vapor,  liquid 

2276 

OUTPUT  709! ‘AR  AF0  ALS  VRS* 

2280! 

Sanpla  each  tharnocoupla  20  tines  and  report 

rnocouplei 

1  vapor-0, l,2i  liquid-3&4 

2284 

FOR  1-0  TO  5 

2288 

Sun-0 

2292 

OUTPUT  709; ‘AS  SA‘ 

2296 

FOR  Ji-1  TO  20 

2300 

ENTER  709iEllq 

2304 

Sun-Sun4'Eliq 

2308 

NEXT  Jl 

2312 

Enf(n-Sun/20 

2316 

T(I)-FNTv8v(Enf(I)) 

2320 

NEXT  I 

2324! 

Conpute  average  tenperature  of  liquid 

2328 

Tlav-(T(3>+T(4))*.5 

2332! 

Conpute  average  tenperature  of  vapor 

2333 

Tvavl-(T(0>+T( 1 >>/2 

2334 

Tvav2-T(2) 

2336 

Tvav-(T(0)+T(l  )+T(2))/3 

2340  . 

IF  ABS(Tlav-Dtld)>.2  THEN 

2344 

IF  Tlav>Dtld  THEN 

2348 

BEEP  4000,. 2 

2352 

ELSE 

2356 

BEEP  250,. 2 

Tld2 


TIbb 
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2360 

2364 

2368 

2372 

2376 

2380 

2384 

2388 

2392 

2396 

2400 

2404 

2408 

2412 

2416 

V 

2420 

2424 

2428 

2432! 

2436! 

2440 

2444 

2448 

2452 

2456 

2460 

2464 

2468 

2472 

2476 

2480 

2484 

2488 

2492 

2496 

2500 

2504 

2508 

2512 

2516 

2520 

2524 

2528 


END  IF 

ELSE 

IF  ABS(Tlav-Dtld».l  THEN 
IF  Atld>Dtld  THEN 
BEEP  3000,. 2 

ELSE 

BEEP  800,. 2 
END  IF 
END  IF 
END  IF 

Errl-Tlav-Oldl 

Oldl-Tlav 

Err2“Tvav-01d2 

01d2^Tvav 

PRINT  USINS  **X.7(MDDD.DD.3X>-iDtld.T(3).T(4>.T(S).Tvavl.Tvav2.Tla 

URIT  2 
60T0  2252 
END  IF 

TAKE  DATA  IF  In-0  LOOP 
IF  Ikol-1  THEN  2452 
BEEP 
Ikol-1 

OUTPUT  709i “AR  AF0  AL5  VRS* 

FOR  1-0  TO  5 

OUTPUT  709 i* AS  SA* 

Sun-0 

FOR  Ji-1  TO  20 
ENTER  709 iE 
Sun-5un+E 

IF  I>2  THEN  EtUi*l>»E 
NEXT  Jl 
Kdl-0 

IF  I >2  THEN 
Eava-Sun/20 
Sun-0. 

FOR  Jk-0  TO  19 

IF  ABS<Et(Jk)-Eeve)<5.eE-6  THEN 
Sun-Sun-t'EK  Jk ) 

ELSE 

Kdl-Kdl-t-l 
END  IF 
NEXT  Jk 


2532 

2536 

2540 

2544 

2548 

peat  data  set’*' 


IF  I>2  THEN  PRINT  USING  ’4X,”Kdi  -  ”,DD"iKdl 
IF  Kdl>10  THEN 
BEEP 
BEEP 

PRINT  USINS  ■4X,"Too  nuch  scattering  in  data  -  re 


2552  goto  1980 


2556 

END  IF 

2560 

END  IF 

2564 

Enf< I  )-Sun/<20-Kdi  ) 

2568 

NEXT  I 

2572  , 

OUTPUT  709i’AR  AF40  AL69  VRS 

2576 

FOR  1-6  TO  35 

2580 

OUTPUT  709 i "AS  SA’ 

2584 

Sun-0 

2588 

FOR  Ji-1  TO  5 
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;S92 
2596 
2600 
2604 
2608 
2612! 
2616! 
2620 
2624 
2628 
2632 
2636 
2640! 
2644! 
2648 
2652 
2656 
2660 
2664 
2668! 
2672 
2676 
2680 
2684 
2688 
2692 
2696 
2700 
2704! 
2708! 
2712 
2716 
2720 
2724 
2728 
2732 
2736 
2740 
2744! 
2748 
2752 
2756 
2760! 
2764 
2768 
2772 
2776! 
ernocoup 
2780 
2784 
2788 
2792 
2796 
2800 
2804 
2808  . 
2812 
2816 
2820 
2821 


ENTER  709 ;E 
Sum«Sun+E 
NEXT  Jl 
Enf < I  )“Sun/5 
NEXT  I 

READ  V0LTA6ES  ( 27-Inat ,28-Sln,29-Aux ) 

OUTPUT  709: ”AR  AF27  AL29  VRS" 

FOR  >0  TO  2 

OUTPUT  709; "AS  SA” 

ENTER  709iUolt<I) 

NEXT  I 

READ  CURRENTS  (30-34-Inst  tubes :35-39-ACTIUE  Dunny ) 

OUTPUT  709: ‘AR  AF30  AL39  VRS' 

FOR  1-0  TO  9 

OUTPUT  709: “AS  SA“ 

ENTER  709:Anp(I> 

NEXT  I 

Read  Currents! 25-Aux  anps, 26-Sin  anps) 

OUTPUT  709: ”AR  AF25  AL26  URS“ 

FOR  1-10  TO  11 

OUTPUT  709; ‘AS  SA* 

ENTER  709;Anp<l) 

NEXT  1 
ELSE 

ENTER  OFile2;Enf<» ),Uolt(*>.Anp(*> 

ENO  IF 

CONVERT  EMF’S  TO  TEMP , UOLT , CURRENT 
FOR  1-0  TO  35 

T<I)-FNTvav(Enf{IJ> 

IF  1>S  ANO  ldtc>0  THEN 
FOR  li-0  TO  ldtc-1 

IF  Ldtcdi  >-1-5+39  THEN  T(l)— 99.99 
NEXT  li 
END  IF 
NEXT  1 

Ntc-nr  of  thernocouples 
Ntc-6 

FOR  li-0  TO  4 

0(  li )-60-Volt<0  )»Anp( li  ) 

Twa-Average  tenperature  of  the  wall 

Twa(li)-0 

Ndte-0 

FOR  1-1  TO  Ntc 

Nn  is  counter  in  tenp  array,  start  at  6  (this  is  the  first  th 
le  in  the  tube  bank  ) 

Nn-I i*6+l+5 
IF  ABS(T(Nn))  99  THEN 
T(Nn  )— 99.99 
Ndtc-Ndtc+1 

ELSE 

Twa( li  )-Twa< li )+T<Nn) 

END  IF 
NEXT  1 

Twa(  li )-Twa( Ii  )/(6-Ndtc  ) 

NEXT  li 

Tlav-(T(3)+T(4))/2 
Tvav«T ( 2 
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Tvav(  T(0  )+T(  1  )+T<2  )  )/3 


2824' 

2828' 

2829! 

2832 

2836 

2840! 

2944! 

2848 

2852 

2856 

2860 

2864 

2868 

2872! 

2876 

2880 

2884! 

2888! 

2892 

2896 

2900 

2904 

2908 

2912 

2916 

2920 

2924 

2928! 

2932! 

2936! 

2940 

2944 

2948 

2952 

2956 

2960 

2964 

2968 

2972 

2976 

2980 

2984 

2988! 

2992! 

2996 

3000 

3004 

3008 

3012 

3016 

3020 

3024 

3028 

3032 

3036 

3040 

3044 

3048 

3052 

3056 


Tlav-T(S) 

Tcu«Tuia(  0  > 

Kcu-FNKcu(Tou)  ! THERMAL  CONDUCTIVITY  OF  COPPER 

!IF  CURVE  FIT  NOT  AVAIL  USE  ARRAY  KCU( •  ) 
FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NE6LECTEO 
FOR  1-0  TO  4 

Tu(  I  )-Twa(  I  )-0(  I  )»LOS(  02/01  )/<  2*PI*Kcu*L  ) 

IF  Ilqv-0  THEN  Texs-Tlav 
IF  IlQv-1  THEN  Taxa-Tvav 
IF  nqv-2  THEN  Tex9-(  Tlav+T(  2  >  )•  .5 
IF  Corr-1  THEN  Thetab( I  )-Ty( I  )-Texs 

IF  Corr-0  THEN  Thetab( 1  )-Tu( I  )-(Texs+.0S6+I».129>  !R-114 

IF  Corr-0  THEN  Thetab( I >-Tw( 1 )-( Tex9+.054+l*. 144 >  !R-113 

NEXT  I 


COMPUTE  VARIOUS  PROPERTIES 


Tfili<i-<Tw<0)+Tex9  )«.S 
Rha-FNRho<  Tf i 1b ) 
Mu-FNMu(TfilB) 
K-FNK<Tf ilB) 
Cp-FNCpdfilm) 
Beta-FNBata( Tf i 1b ) 
Ni-Mu/Rho 
Alpha-K/(Rho*Cp ) 
Pr-Ni/Alpha 


!FILM  TEMPERATURE 
•DENSITY 
! VISCOSITY 

! THERMAL  CONDUCTIVITY 
! SPECIFIC  HEAT 
! THERMAL  EXPANSION 
!KINEMATIC  VISCOSITY 
•THERMAL  DIFFUSIVITV 
IPRANDTL 


COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 
FOR  UNENHANCEO  ENO<S) 

Lu-Lua(Itt) 

Hbar-190 

Fa-<Hbar*P/<Kcu*A)  >“.S*Lu 
Tanh-FNT8nh(Fa ) 

Theta<  0  )-Thatab<  0  )»Tanh/Fe 

Xx*< 9.81*Bete*Thetab<0  )»0o'‘3*Tanh/( Fa*Ni*Alpha  >  >* . 166667 
Yy-<l+( .5S9/Pr)'<9/16))*(8/27> 

Hbarc-K/Do*< .6+.387»Xx/Yy  )'2 
IF  ABS<<Hbar-Hbarc)/Hbarc».001  THEN 
Hbar-<Hbar+Hbarc  )*.5 
SOTO  2948 
ENO  IF 


COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCEO  END(S) 
Ql<  0  )-<  Thetab<  0  )*Tanh  )•( ( Hbar*P»Kcu»A  )" .5  ) 
Oq-Ol<0>+Oq 
Z-Z+1 

IF  Z-1  THEN 
Lu-Lv 
GOTO  2944 
ENO  IF 
Z-0 

01pct-0q/0( 0 ) 

Oq-0 

A9-PI*02»L 
,  FOR  Ii-0  TO  4 

Ql(Ii)-Qlpct»Q(Ii > 

Odpdl  )-<0(Ii  >-Ql(Ii)>/A9 
Htube( Ii  )-Odp< Ii >/Thetab( li  ) 

NEXT  Ii 
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PRINTER  IS  701 


3060 
30641 

3068'  RECORD  TIME  OF  ORTA  TAKING 

3072  IF  In-0  THEN 

3076  OUTPUT  709i"TO“ 

3080  ENTER  709iTold$ 

3084  END  IF 

30881 

3092!  OUTPUT  DATA  TO  PRINTER 

3096  PRINTER  IS  701 

3100  PRINT 

3104  PRINT  USING  "10X,‘"Oata  Set  Nunber  -  * ' ,DDO,2X , UA“ : J .Told* 

3108  PRINT 

3112  PRINT  USING  "10X.”  Tvl  Tv2  Tv3  Tldl  Tld2  Tld3 

Tvav  Tldav 

3116  PRINT  USING  "IGX .8< MDD.DD,2X >" ; T( 0  ),T< 1 >,T( 2 >.T( 3  )  ,T< 4  ).T( S > ,Tvav . Tlav 

3120  PRINT 

3124  PRINT  USING  ‘EX.^Tube  Uall  Tenperatures  (Deg  C>  Tnave  (}dp 

H  Thetab""" 

3128  PRINT  USING  "6X,'‘*  1  2  3  4  5  G  (Deg  C>  (U/n* 

2)  <W/n*2.K)  (K)“““ 

3132  Jj-0 

3136  FOR  Ii-0  TO  Nht-1 

3140  FOR  Ji-0  TO  S 

3144  Tp( Ji >-T( Ii»S+Jj+6 > 

3148  Jj-Jj+1 

3152  NEXT  Jx 

3156  Jj-Ii+1 

31S0  FOR  Jx-0  TO  4 

31G4  Tn(Ji)-l+Ji 

3168  NEXT  Ji 

3172  PRINT  USING  “SX ,D. 1X,7(MDO.OD >, 1X,2( MZ.3DE > , IX , KMDO.DD >* i Tn(  Ii ) ,T 

p<0),Tp(l ),Tp<2>,Tp(3),Tp(4),Tp(S),Tue<Ii>,0dp(Ii >,Htube(Ii ),Thetab(Ii > 

3176  NEXT  li 

3180  Ok-1 

3184  IF  In-0  THEN 

3188  BEEP 

3192  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( 1-Y( default >,0-N >?■ ,0k 

3196  END  IF 

3200!  J-the  counter  for  date  seta 

3204  IF  Ok-1  OR  In-1  THEN  J-J+1 

3208  IF  Ok-1  AND  In-0  THEN  OUTPUT  *Filel ;Enf( * >,Uolt( *  ),Anp( •  ) 

3212  IF  In-1  OR  Ok-1  THEN  OUTPUT  #Plot iOdp( •  ),Htube( •  ) , Thetab( • ) 

3216  Go_on-1 

3220  IF  In-0  THEN 

3224  BEEP 

3228  INPUT  "WILL  THERE  BE  ANOTHER  RUN  ( 1-Y' default  ).0-N>’" ,Go_on 

3232  Nrun-J 

3236  IF  Go_on-0  THEN  3272 

3240  IF  So_on' '0  THEN  Repeat 

3244  ELSE 

3248  IF  J  Nrun+1  THEN  Repeat 

3252  END  IF 

3256  St-1 

3260  BEEP 

3264  .  INPUT  "ARE  YOU  SURE  YOUR  READY  TO  TERMINATE  ( 1-Y( DEFAULT », 0-NO St 

3268  Go_on-l 

3272  IF  St>0  THEN  3280 

3276  IF  St-0  THEN  GOTO  3240 

3280  IF  In-0  THEN 
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3284  BEEP 

3288  PRINT 

3292  PRINT  USING  "lOX. ’"NOTE-  ‘“.ZZ,"  data  runs  uere  stored  in  file  ” 

,10A"; J-l.D2file$ 

3296  ASSIGN  SFilel  TO  • 

3300  OUTPUT  0File2:Nrun-l 

3304  ASSIGN  0Filel  TO  OlfileS 

3308  ENTER  0Fi lei : DateS ,Ldtc< *  )| Itt ,6op ,Nht ,Natp ,Nrt ,Corr 

3312  OUTPUT  0File2iOateS,Ldtc( * >, It t ,Bop ,Nht ,Natp ,Nrt ,Corr 

331G  FOR  I-l  TO  Nrun-l 

3320  ENTER  @Filel :Enf< *  ),Uolt( •  ),Anp( *  ) 

3324  OUTPUT  0Fi le2 jEnf ( •  ) .Uol t< • >, AnpC •  ) 

3328  NEXT  I 

3332  ASSIGN  0Filel  TO  • 

333G  PURGE  •DUMMY" 

3340  END  IF 

3344  BEEP 

3348  PRINT 

3352  PRINT  USING  “10X,""NOTE:  '“.ZZ,"  X-Y  pairs  yere  stored  in  plot  data  f 

ile  ".10A“!j-l,Pfile$ 

3356  ASSIGN  0File2  TO  • 

3360  ASSIGN  0Plot  TO  • 

3364  BEEP 

3368  SUBENO 
3372! 

337G!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

3380! 

3384  DEF  FNKcu(Tcu> 

3388!  OFHC  COPPER 

3392  Tk-Tcu+273.15  !C  TO  K 

3396!  Kcu-434-.112»Tk  !250-300K  USE  FOR  R-114  #2.2  C 

3400  Kcu>433.0-.l*Tk  !200-400K  USE  FOR  R-113  047.5  C 

3404  RETURN  Kcu 
3408  FNENO 
3412! 

3416  DEF  FNMu(T) 

3420!  CURVE  FIT  OF  VISCOSITY 
3424  Tk-T+273.15  !C  TO  K 

3428!  Mu-EXP(-4.4636+(1011.47/Tk))*1.0E-3  IR-114  170-360  K 

3432  Mu-.0000134*<10“(503/(Tk-2.15)>)  !R113 

3436  RETURN  Mu 

3440  FNENO 

3444! 

3448  DEF  FNCp<T) 

3452!  CURVE  FIT  OF  Cp 

3456  Tk-Tf273.15  !C  TO  K 

34601  Cp-. 40188+1. 65007E-3»Tk+1.51494E-6»Tk'2-6.G7853E-10»Tk"3  iR-114  180-400  K 

3464  Cp-( 929+1. 03»T)». 001  'R-113 

3468  Cp-Cp*1000 

3472  RETURN  Cp 

3476  FNENO 

3480! 

3484  DEF  FNRho(T) 

3488  Tk-T+273.15  'C  TO  K 

3492  X-l-a.8»Tk/7S3.95)  'K  TO  R 

3496'  Rp-36. 32+61. 146414«X-( 1/3  )+16 .418015-X+17. 476838»X ' .5+1 . 1 19828«X'2 
35001  Ro-Ro/ .062428  iR-114 

3504  Ro-1.6207479E+3-T»(:.2186346+T»2.3578291E-3)  iR-113 

3508  RETURN  Ro 
3512  FNENO 


159 


3516' 

35:B  DEF  FNPra)  'GOOD  FOR  R-lU/R-113 
3524  Pp-FNCp(T)»FNf1u(T)/FNK(T) 

3528  RETURN  Pr 
3532  FNEND 
35361 

3540  DEF  FNK(T) 

35441  T>'360  K  WITH  T  IN  C 
3548  K-.071-.000261»T 

3552  RETURN  K 
3556  FNEND 
3560! 

3564  DEF  FNTanh<Fe> 

3568  P-EXP(Fe) 

3572  Q-EXP(-Fe) 

3576  Tanh-(P-0)/(P+Q) 

3580  RETURN  Tanh 
3584  FNEND 
35881 

3592  DEF  FNTvsvCU) 

3596  COM  /Cc/  C(7) 

3600  T-C<0> 

3604  FOR  >1  TO  7 
3608  T-T+C(I  I'U-I 
3612  NEXT  I 
3616  RETURN  T 
3620  FNEND 
3624! 

3628  DEF  FNBeta(T> 

3632  Rop-FNRho<T+.l ) 

3636  Ron-FNRho<T-.l ) 

3640  Beta“-2/<Rop+Ron )»(Rop-RoB)/.2 

3644  RETURN  Bata 

3648  FNEND 

3652  DEF  FNPoly(X) 

3656  COM  /Cply/  A< 10. 10 >.C< 10 ).B<4 ),Nop, Iprnt .Opo, Ilog 
3660  Xl-X 
3664  Poly-B<0) 

3668  FOR  I-l  TO  Nop 

3672  IF  Ilog-1  THEN  X1-L06(X) 

3676  Poly-Poly+B<I  )»X1“I 

3680  NEXT  I 

3684  IF  Ilog-1  THEN  Poly-EXP< Poly  ) 

3686  RETURN  Poly 
3692  FNEND 
36961 

37001XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

3704! 

3708  SUB  Poly(Ofile$C  ),Np,Itn) 

3712  DIM  R(10>,S<10),Sy<12!,S<(12),Xr.<l00>.Yy(100),Xy(17) 

3716  COM  /Cply/  A( 10, 10  )  .C( 10  )  ,B( 4  )  ,N, Iprnt ,0po. I  log 
3720  COM  /Xxyy/  Xpf5),Yp(5) 

3724  FOR  1-0  TO  4 
3728  B(  I  >-0 
3732  NEXT  I 
3736  Iyi-1 
3740  BEEP 

3744  INPUT  "ENTER  DATA  FILE  NAME" ,Df i leS(  0 ) 

3748  BEEP 

3752  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS", Np 
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3756  BEEP 

3760  INPUT  "LIKE  TO  EXCLUDE  DATA  PAIRS  ( 1»Y,0«N(DEFAULT ; )’•  led 
3764  IF  Ied-1  THEN 
3768  BEEP 

3772  INPUT  "ENTER  NUMBER  OF  PAIRS  TO  BE  EXCLUDED" , Ipex 
3776  END  IF 

3780  ASSISN  •File  TO  Ofilet(0) 

3784  N-2 

3788  BEEP 

3792  INPUT  “ENTER  THE  ORDER  OF  POLYNOMIAL  (DEFAULT-2)  ”  N 

3796  FOR  1-0  TO  N 

380P  Sy<I)-0 

3804  Sx(I)-0 

3808  NEXT  I 

3812  IF  Ied-1  AND  In-l  THEN 
3816  FOR  1-1  TO  Ipex 

3820  ENTER  •File;Xy<«) 

3824  NEXT  I 

3828  END  IF 

3832  FOR  r-1  TO  Np-Ipex 

3836  ENTER  •File:Xy<*) 

3840  IF  Opo-0  THEN 

3844  Y-Xy(Itn-l) 

3848  X-Xy(ll+Itn) 

3852  END  IF 

3856  IF  Opo-1  THEN 

3860  Y-Xy<5+Itn) 

3864  X-Xy(ll+Itn) 

3868  END  IF 

3872  IF  Opo-2  THEN 

3876  Y-Xy<S+Itn) 

3880  X-Xydtn-l) 

3884  END  IF 

3888  IF  1  log-1  THEN 

3892  X-L06(X> 

3896  Y-L06(Y> 

3900  END  IF 

3904  Xx<I)-X 

3908  Yy(  I  )-Y 

3912  R<0)-Y 

3916  Sy(0)-Sy(0>'»'Y 

3920  S(  1  )-X 

3924  Sx(l  )-Sx(l»+X 

3928  FOR  J-1  TO  N 

3932  R<  J  )-R<  J-1  )»X 

3936  Sy(J)-Sy(J)+R(J) 

3940  NEXT  J 

3944  FOR  J-2  TO  N<*2 

3948  S<  J  )-S(  J-1  )»X 

3952  S/.(J)-Sx(J>+S(J) 

3956  NEXT  J 

3960  NEXT  I 
3964  Sx(0)-Np 
3968  FOR  1-0  TO  N 
3972  C(  I  )-Sy<  I ) 

3976  ,  FOR  J-0  TO  N 

3980  A< I.J  )-Sx< I+J  ) 

3984  NEXT  J 

3988  NEXT  I 

3992  FOR  1-0  TO  N-1 
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3996  CAUDivide(I) 

4000  CALL  Subtractf  I  +  n 

4004  NEXT  I 

4008  B(N)>C(N)/A(N,N) 

4012  FOR  1-0  TO  N-1 
4016  )-C<N-l-I ) 

4020  FOR  J-0  TO  I 

4024  B(N-1-I  )-B<N-l-I  )-A( N-l-I ,N-J  )«B(N-J  > 

4028  NEXT  J 

4032  B<N-1-I  )-B(N-l-I  )/A( N-l-I .N-l-I ) 

4036  NEXT  I 

4040  IF  Iprnt-0  THEN 

4044  PRINT  USING  " 12X 'EXPONENT  COEFFICIENT'" 

4048  FOR  1-0  TO  N 

4052  PRINT  USING  ' ISX ,D0,5X ,M0.7DE' : I ,B( I > 

4056  NEXT  I 

4060  PRINT  “  ' 

4064  PRINT  USING  '12X,"0ATA  POINT  X  Y  Y<  CALCULATED )  01 

SCREPANCY'" 

4068  FOR  I-l  TO  Np 

4072  Yc-B(0) 

4076  FOR  J-1  TO  N 

4080  Yc-Yc+B( J  )*Xx< I >'J 

4084  NEXT  J 

4088  0-Yy< I  )-Yc 

4092  PRINT  USING  ' 1SX.30,4X.4(MD.SDE, IX >' : I .Xx( I > , Yy( I ), Yc ,D 

4096  NEXT  I 

4100  END  IF 

4104  ASSIGN  0File  TO  • 

4108  SUBENO 
4112) 

4116  SUB  Oivide(M) 

4120  COM  /Cply/  A( 10, 10 ),C( 10),B(4 ),N, Iprnt ,0po. Ilog 
4124  FOR  I-M  TO  N 
4128  Ao-A(I,n) 

4132  FOR  J-fl  TO  N 

4136  A(I,J>-A<I,J)/Ao 

4140  NEXT  J 

4144  C(I)-C(I)/Ao 

4148  NEXT  1 

4152  SUBENO 

4156) 

4160  SUB  SubtracKK) 

4164  COM  /Cply/  A(10. 10), C< 10), B<4),N, Iprnt. Opo, Ilog 

4168  FOR  I-K  TO  N 

4172  FOR  J-K-1  TO  N 

4176  A(I,J)-A(K-1,J)-A(I,J) 

4180  NEXT  J 

4184  C<  I  )-C(K-l  >-C(  I  ) 

4188  NEXT  I 
4192  SUBENO 
4196) 

4200  SUB  PI  in 

4204  COM  /Cply/  A( 10, 10  ) , C< 10  )  .B( 4  ) ,N. Iprnt , Opo, I  log 
4208  COM  /Xxyy/  Xx(5>,Yy(S) 

4212  PPINTER  IS  70S 
4216  BEEP 

4220  INPUT  "SELECT  (0-h/h0X  sane  tube, l-h( HF >/h( sn >' , Irt 
4224  BEEP 

4228  INPUT  "WHICH  Tsat  ( l-B.7,0— 2.2  )',Isat 
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4232  Xi<iin»0 
4236  Xna.x-ia 
4240  Xstep-2 
4244  IF  Irt-0  then 
4248  Ynin“0 

4252  Ypiax-1.4 

4256  Y3tep-.2 

4260  ELSE 

4264  Ynin“0 

4268  Yi>iai<*l5 

4272  Ystep-5 

4276  END  IF 
4280  BEEP 

4284  PRINT  "INiSPlilP  2300,2200.8300.6800!” 

4288  PRINT  “SC  0, 100,0. 100 i TL  2,0:“ 

4292  Sfx«100/< Xnax-Xnin ) 

4296  Sfy“100/(  Ynax-Ytnin ) 

4300  PRINT  “PU  0,0  PD“ 

4304  FOR  Xa-Xnin  TO  Xnax  STEP  Xstep 

4308  X-(Xa-Xiiin  )»Sfx 

4312  PRINT  “PA"iX,“,0:  XT:  “ 

4316  NEXT  Xa 

4320  PRINT  "PA  100,0:PUi“ 

4324  PRINT  “PU  PA  0,0  PD“ 

4328  FOR  Ya-Ynin  TO  Yoiax  STEP  Ystep 
4332  Y-( Ya-Ynin  )»Sfy 

4336  PRINT  “PA  0,*iY,“YT* 

4340  NEXT  Ya 

4344  PRINT  "PA  0,100  TL  0  2“ 

4348  FOR  Xa”Xnin  TO  Xnax  STEP  Xstep 

4352  X-(Xa-Xnin)»Sfx 

4356  PRINT  “PA“ i X , “ , 100;  XT” 

4360  NEXT  Xa 

4364  PRINT  “PA  100,100  PU  PA  100,0  PD“ 

4368  FOR  Ya«Ynin  TO  Ynax  STEP  Ystep 

4372  Y-( Ya-Ynin  )*Sfy 

4376  PRINT  “PD  PA  100,“,Y.-YT" 

4380  NEXT  Ya 

4384  PRINT  “PA  100,100  PU“ 

4388  PRINT  "PA  0,-2  SR  1.5,2“ 

4392  FOR  Xa"Xnin  TO  Xnax  STEP  Xstep 
4396  X-<Xa-Xnin)*Sfx 

4400  PRINT  “PA’;X.“,0:“ 

4404  PRINT  “CP  -2,-l;LB”;Xa: ■“ 

4408  NEXT  Xa 

4412  PRINT  "PU  PA  0,0“ 

4416  FOR  Ya»Ynin  TO  Ynax  STEP  Ystep 
4420  IF  ABS(Ya)<l.E-5  THEN  Ya-0 

4424  Y»( Ya-Ynin  )»Sfy 

4428  PRINT  "PA  0,"!Y,"“ 

4432  PRINT  "CP  -4,-.25.LB";Yai "■ 

4436  NEXT  Ya 

4440  Xlabelf“"Oil  Percent" 

4444  IF  Irt-0  THEN 
4448  YlabelS--h/h0S“ 

4452  -  ELSE 

4456  Ylabel*”"h/hsnooth" 

4460  END  IF 

4464  PRINT  "SR  1.5.2:PU  PA  50,-10  CP" ; -LEN( X labels '/2 ; "0 : LB X labe !$ 
4468  PRINT  "PA  11,50  CP  0.";-LEN< yiabel$)/2»5/6. "01  0 , 1 . LB " ^ Y1 abe 1  $ 
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4472  PRINT  “CP  0,0“ 

4476  Ipn-0 
4480  BEEP 

4484  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  ( 1-Y , 0-N , Okp 

4488  Icn-0 

4492  IF  Okp-1  THEN 

4496  BEEP 

4500  INPUT  “ENTER  THE  NAME  OF  THE  DATA  FILE",D  fileS 
4504  BEEP 

4508  INPUT  “SELECT  (0-LlNEAR.  I-LOG( X , Y )" , I  log 
4512  ASSIGN  eFile  TO  0  fileS 
4516  BEEP 

4520  INPUT  “ENTER  THE  BEGINNING  RUN  NUMBER”, Md 
4524  BEEP 

4528  INPUT  “ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" , Npairs 
4532  BEEP 

4536  INPUT  “ENTER  DESIRED  HEAT  FLUX".0 

4540  BEEP 

4544  PRINTER  IS  1 

4548  PRINT  USING  “4X,“ “Select  a  symbol :“““ 

4552  PRINT  USING  "4X.““1  Star  2  Plus  sign"““ 

4556  PRINT  USIN6  "4X,“*3  Circle  4  Square" ““ 

4560  PRINT  USING  "4X,“"S  Ronbus”*” 

4564  PRINT  USING  “4X,““G  Right-side-up  triangle""" 

4568  PRINT  USING  *4X,““7  Up-side-douin  triangle" ““ 

4572  INPUT  Sym 

4576  PRINTER  IS  70S 
4580  PRINT  "PU  01" 

4584  IF  Syrt-1  THEN  PRINT  “SM*” 

4588  IF  Syit-2  THEN  PRINT  “SM+“ 

4592  IF  Syn-3  THEN  PRINT  “SMo“ 

4596  Nn*4 

4600  IF  Ilog-l  THEN  Nn-1 
4604  IF  Md>l  THEN 
4608  FOR  I-l  TO  (Md-1) 

4612  ENTER  0File:Xa,Ya 

4616  NEXT  I 

4620  END  IF 
4624  Ql-0 

4S28  IF  Ilog-l  THEN  O-LOG(Q) 

4632  FOR  I-l  TO  Npairs 
4636  ENTER  9FileiXa,B( • > 

4640  Ya-B(0) 

4644  FOR  K-1  TO  Nn 

4648  Ya-Ya+B<K  )«0“K 

4652  NEXT  K 

4656  IF  Ilog-l  THEN  Ya-EXP(Ya) 

4660  IF  IIOO-0  THEN  Ya-01/Ya 

4664  IF  Irt-0  THEN 

4668  IF  Xa-0  THEN 

4672  Yo-Ya 

4676  Ya-1 

4680  ELSE 

4684  Ya-Ya/Yo 

4688  END  IF 

4692  ■  ELSE 

4696  Hsn-FNHsnooth<0,Xa, Isat > 

4700  Ya-Ya/Hsn 

4704  END  IF 

4708  I-l  )-Xa 
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4712  yy(I-l)-Ya 

471B  X»( Xa-Xnin  )*Sf K 

4720  Y-<Ya-Ynln>»Sfy 

4724  IF  SyM>3  THEN  PRINT  ‘SM‘ 

4728  IF  Sy«<4  THEN  PRINT  "SR  1.4, 2. 4" 

4732  PRINT  "PA'.X.Y,"" 

4736  IF  SyM>3  THEN  PRINT  ‘SR  1.2, 1.6" 

4740  IF  Sy)0-4  THEN  PRINT  "UC2, 4 . 99, 0, -8 , -4 , 0 ,0 , 8 , 4 , 0,  ;  " 
4744  IF  Syn-S  THEN  PRINT  “UC3.0 , 99, -3 ,-6 .-3 ,6 , 3 ,6 , 3 , -6 : " 
4748  IF  Syoi-e  THEN  PRINT  “UC0,S.  3, 99, 3 , -8 . -6 . 0, 3 . 8 ;  " 

4752  IF  SyM-7  THEN  PRINT  "UC0, -S .3, 99, -3 , 8 , B ,0, -3, -8: “ 
4756  NEXT  I 
4760  BEEP 

4764  ASSIGN  «Fila  TO  • 

4768  END  IF 
4772  PRINT  "PU  SM" 

4776  BEEP 

4780  INPUT  “UANT  TO  PLOT  A  POLYNOMIAL  ( 1-Y ,0-N >? " ,Okp 
4784  IF  Okp-1  THEN 
4788  BEEP 

4792  INPUT  "SELECT  (0-LINEAR, 1-L06(X , Y  ) )" , I  log 

4796  lprnt-1 

4800  CALL  Polydtn) 

4804  FOR  Xa-Xnin  TO  Xnax  STEP  Xstep/25 

4808  Icn”Icn+i 

4812  Ya-FNPoly(Xa) 

4816  Y-<  Ya-Yinin)»Sfy 

4820  X-( Xa-Xnin  )»Sfx 

4824  IF  Y<0  THEN  Y-0 

4828  IF  Y>100  THEN  GOTO  4868 

4832  Pu-0 

4836  IF  Ipn-1  THEN  Idf-Icn  MOD  2 

4840  IF  Ipn-2  THEN  Idf-Icn  MOO  4 

4844  IF  Ipn-3  THEN  Idf-Icn  MOO  8 

4848  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

4852  IF  Ipn-S  THEN  Idf-Icn  MOD  32 

4856  IF  Idf-1  THEN  Pu-l 

4860  IF  Pu-0  THEN  PRINT  "PA" ,X,Y, "PD" 

4864  IF  Pu-l  THEN  PRINT  ■PA",X,Y, "PU" 

4868  NEXT  Xa 

4872  PRINT  "PU* 

4876  Ipn-Ipn+1 

4880  GOTO  4480 

4884  END  IF 
4888  BEEP 

4892  INPUT  "UANT  TO  QUIT  ( 1-Y,0-N>‘?"  ,Iquit 
4896  IF  Iquit-1  THEN  4904 
4900  GOTO  4480 
4904  PRINT  "PU  SP0" 

4908  SUBEND 
4912  SUB  Stats 
4916  PRINTER  IS  701 
4920  J-0 

4924  K-0 

4928  BEEP 

4932  INPUT  "PLOT  FILE  TO  ANALYZE’" ,Fi leS 
4936  ASSIGN  0File  TO  FileS 
4940  BEEP 

4944  INPUT  "LAST  RUN  No’< 0-OUIT  )* ,Nn 
4948  IF  Nn-0  THEN  5092 


165 


495Z  Nn-Nn-j 

4956  Sx-0 

4960  Sy-0 

4964  Sz-0 

4968  Sx5“0 

4972  Sy9-0 

4976  929-0 

4980  FOR  I-l  TO  Nn 

4984  J-J+1 

4988  ENTER  SFileiO.T 

4992  H-Q/T 

4996  Sx-Sx+Q 

5000  Sx9-Sx9-t-0'2 

5004  Sy-Sy+T 

5008  Sy9-Sy9+T“2 

5012  Sz-Sz+H 

5016  Sz9-Sz9+H“2 

5020  NEXT  I 

5024  Qave-Sx/Nn 

5028  Tave-Sy/Nn 

5032  Have-Sz/Nn 

5036  Sdevq-SOR(ABS( (Nn»Sx9-Sx“2  )/(Nn*<Nn-l ) ) > ) 

5040  Sdevt-SQR( ABS( <Nn*Sy9-Sy*2 )/<Nn*(Nn-l ) )>  > 

5044  Sdevh-SQR( ABS< (Nn»Sz9-Sz*2  >/(Nn*<Nn-l ) ) ) ) 

5048  SN-100«Sdevh/Have  . 

5052  Sq-100«Sdevq/Qava 
5056  St-100*Sdevt/Tave 
5060  IF  K-1  THEN  5084 
5064  PRINT 

5068  PRINT  USIN6  "1 IX, ““DATA  FILE" ,14A’ i Filet 
5072  PRINT 

5076  PRINT  USIN6  *1 IX, ’"RUN  Htube  SdevH  Odp  SdevQ  Thetab  SdovT” 

m 

5080  K-1 

5084  PRINT  USIN6  ’ 1 IX ,00, 2(2X,D.30E, IX ,30.20 >,2X ,00.30, IX ,30.20’ ; J, Have, Sh, Gave 

,Sq,Tave,St 

5088  GOTO  4940 

5092  ASSIGN  OFilel  TO  • 

5096  PRINTER  IS  1 
5100  SUBENO 
5104  SUB  Coef 

5108  COM  /Cply/  A( 10, 10 > ,C( 10  )  ,B< 4  ),N. Iprnt ,0po, I  log 
5112  BEEP 

5116  INPUT  "GIVE  A  NAME  FOR  CROSS-PLOT  FILE’.Cpft 
5120  BEEP 

5124  INPUT  "OUTPUT  TYPE  <0-q  va  Dt ,  1-h  v9  Ot,  2-h  v9  q)’,0po 
5128  CREATE  BOAT  Cpff,G 
5132  ASSIGN  »File  TO  Cpft 
5136  BEEP 

5140  INPUT  "SELECT  (0-LINEAR.l-LO6<X,Y> )",Ilog 
5144  BEEP 

5148  INPUT  "ENTER  OIL  PERCENT  < -l-STOP  )' ,Bop 
5152  BEEP 

5156  INPUT  "ENTER  TUBE  NUMBER  <1,  2,  3,  4,  OR  S'-.Itn 
5160  IF  Bop'.0  THEN  5176 
5164  CALL  PolyCItn) 

5169  OUTPUT  0FileiBop,B( * ) 

5172  GOTO  5144 

5176  ASSIGN  0File  TO  • 

5180  SUBENO 
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51841 

5138IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

5192! 

5196  SUB  Plot 

5200  con  /Cply/  A(10,10),C(10),8<4>,Nop.Iprnt,Opo,Ilog 
5204  DIM  Xy(17) 

5208  INTEGER  li 
5212  PRINTER  IS  1 
5216  BEEP 
5220  Idv-1 

5224  input  “LIKE  DEFAULT  VALUES  FOR  PLOT  ( 1-Y( DEFAULT  )  ,0»N  )?", Idv 
5228  Opo-0 
5232  BEEP 

5236  PRINT  USING  “4X. ““Select  Option:“““ 

5240  PRINT  USING  “GX,““0  q  versus  del ta-T( DEFAULT >““ “ 

5244  PRINT  USING  “6X,““1  h  versus  delta-T“““ 

5248  PRINT  USING  “6X,““2  h  versus  q“““ 

5252  INPUT  Opo 
5256  BEEP 

5260  INPUT  “SELECT  UNITS  ( 0-SI( DEFAULT), 1-ENGLISH )“, lun 
5264  PRINTER  IS  705 
5268  IF  IdvOl  THEN 
5272  BEEP 

5276  INPUT  “ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS“,Cx 
5280  BEEP 

5284  input  “ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS“,Cy 
5288  BEEP 

5292  INPUT  “ENTER  MIN  X-UALUE  (MULTIPLE  OF  10)“,Xnxn 
5296  BEEP 

5300  INPUT  “ENTER  MIN  Y-UALUE  (MULTIPLE  OF  10>“,Ynin 
5304  ELSE 

5308  IF  Opo-0  THEN 
5312  Cy-2 

5316  Cx-2 

5320  Xnln-1 

5324  Ynin-1000 

5328  END  IF 

5332  IF  Opo-1  THEN 

5336  Cy-2 

5340  Cx-2 

5344  Xnin-1 

5348  Ynin-100 

5352  END  IF 

5356  IF  Opo-2  THEN 

5360  Cy-2 

5364  Cx-2 

5368  Xnin-1000 

5372  Ymin-100 

5376  END  IF 

5380  END  IF 
5384  BEEP 

5388  PRINT  "IN;SPliIP  2300,2200.8300,6800;“ 

5392  PRINT  "SC  0, 100,0, 100 ; TL  2,0;" 

5396  Sfx-100/Cx 
5400  Sfy-100/Cy 
5404  BEEP 

5408  INPUT  “WANT  TO  BY-PASS  CAGE  <1-Y,  0-NC( DEFAULT )", Ibp 
5412  IF  Ibp-I  THEN  5908 
5416  PRINT  "PU  0,0  PD“ 

5420  Nn-9 
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5424  FOR  I-l  TO  C.^+i 
5408  Xat«Xi>ii,n»10*(  I-l  ) 

5432  IF  I-Cx+1  THEN  Nn-1 

5436  FOR  J-1  TO  Nn 

5440  IF  J-1  THEN  PRINT  "TL  2  0’ 

5444  IF  J-2  THEN  PRINT  'TL  1  0" 

5448  Xa-Xat»J 

5452  X-L6T<Xa/Xi>iin)«Sfx 

5456  PRINT  “PA';X,“,0i  XT:“ 

5460  NEXT  J 

5464  NEXT  I 

5468  PRINT  “PA  100,0:PU:' 

5472  PRINT  ”PU  PA  0,0  PO“ 

5476  Nn-9 

5480  FOR  I-l  TO  Cy+1 
5484  Yat-Yi>iin»10“<  I-l  ) 

5488  IF  I-Cy+1  THEN  Nn-1 

5492  FOR  J-1  TO  Nn 

5496  IF  J-1  THEN  PRINT  "TL  2  0‘ 

5500  IF  J-2  THEN  PRINT  "TL  1  0“ 

5504  Ya-Yat«J 

5508  Y-LGT<Ya/YKiin>»Sfy 

5512  PRINT  *PA  0,"!Y,*YT" 

5516  NEXT  J 

5520  NEXT  I 

5524  PRINT  ‘PA  0,100  TL  0  2" 

5528  Nn-9 

5532  FOR  I-l  TO  Cx+1 
5536  Xat-Xnin»10"<I-l ) 

5540  IF  I-Cx+1  THEN  Nn-1 

5544  FOR  J-1  TO  Nn 

5548  IF  J-1  THEN  PRINT  "TL  0  2" 

5552  IF  J>1  THEN  PRINT  "TL  0  1" 

5556  Xa-Xat*J 

5560  X-LGT<Xa/Xinin)*Sfx 

5564  PRINT  "PA" iX, * ,100;  XT" 

5568  NEXT  J 

5572  NEXT  I 

5576  PRINT  "PA  100,100  PU  PA  100,0  PD" 

5580  Nn-9 

5584  FOR  I-l  TO  Cy+1 

5588  Yat-Ynin»10“( I-l > 

5592  IF  I -Cy+1  THEN  Nn-1 

5596  FOR  J-1  TO  Nn 

5600  IF  J-1  THEN  PRINT  "TL  0  2" 

5604  IF  J>1  THEN  PRINT  "TL  0  1" 

5608  Ya-Yat»J 

5612  Y-L6T(Ya/Ynin>»Sfy 

5616  PRINT  ”PD  PA  100,",Y,’YT- 

5620  NEXT  J 

5624  NEXT  I 

5628  PRINT  “PA  100,100  PU" 

5632  PRINT  "PA  0,-2  SR  1.5,2" 

5636  Ii-LGT(xnin) 

5640  FOR  I-l  TO  Cx+1 

5644  *  Xa-Xt>iln«10"(  I-l  ) 

5648  X-LGT( Xa/Xnin >oSfx 

5652  PRINT  "PA-iX.'^O. " 

5656  IF  Ii,’-0  THEN  PRINT  "CP  -2,-2:LB10;PR  - 

5660  IF  li.0  THEN  PRINT  "CP  -2.-2;LB10:PR  0, 
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5664  Ii-Ii+1 

5668  NEXT  I 

5672  PRINT  *PU  PA  0,0‘ 

5676  Ii-L6T(Ymin) 

5680  Y10-10 

5684  FOR  I-l  TO  Cy+1 

5688  Ya-Ynin»10*< I-l ) 

5692  Y-L6T(Ya/YMin>»Sfy 

5696  PRINT  “PA  0,’;Y,““ 

5700  PRINT  “CP  -4,-.25:LB10iPR  -2,2:LB“ : li ; ” “ 

5704  Ii”Ii+l 

5708  NEXT  I 
5712  BEEP 
5716  Idl-1 

5720  INPUT  “WANT  USE  DEFAULT  LABELS  ( 1-Y(  DEFAULT  >,  0-N  Idl 
5724  IF  IdlOl  THEN 
5728  BEEP 

5732  INPUT  “ENTER  X-LABEL” .Xlabel* 

5736  BEEP 

5740  INPUT  “ENTER  Y-LABEL" , Ylabel* 

5744  END  IF 

5748  IF  0po<2  THEN 

5752  PRINT  “SR  1,2:PU  PA  40,-14;“ 

5756  PRINT  “LBCTjPR  -1,6,3  PO  PR  1,2,0  PU:PR  .S,-4:LBuo;PR  .5,1:“ 

5760  PRINT  “LB-T;PR  .5. -1 :LBsat ,PR  .5,1;" 

5764  IF  Iun-0  THEN 

5768  PRINT  "LB)  /  <K)" 

5772  ELSE 

5776  PRINT  "LB)  /  <F)" 

5780  END  IF 

5784  END  IF 
5788  IF  Opp-2  THEN 
5792  IF  lun-0  THEN 

5796  PRINT  "SR  1.5,2:PU  PA  40,-14:LBq  /  (U/m:SR  1,1.5:PR  0.S,l!LB2:SR  1 

.S,2!PR  0.5,-1;LB)" 

5800  ELSE 

5804  PRINT  “SR  1.5,2;PU  PA  34,-14:LBq  /  (Btu/hr;PR  .5,.5:LB.:PR  .5, -.5; 

■ 

5808  PRINT  "L8ft;PR  .5,1:SR  l,l.S:LB2iSR  1.5,2iPR  .5,-l;LB):" 

5812  END  IF 

5816  END  IF 
5820  IF  Opo-0  THEN 
5824  IF  Iun-0  THEN 

5828  PRINT  "SR  1.5,2:PU  PA  -12,40;OI  0,l:LBq  /  (W/in:PR  -1,0.5:SR  1,1. 5;L 

B2;SR  1.5,2:PR  1,.5;LB)“ 

5832  ELSE 

5836  PRINT  "SR  1.5,2:PU  PA  -12,3:;0I  0,l;LBq  /  (Btu/hr;PR  -.5, .5;LB. ;PR 

.5..5:” 

5840  PRINT  "LBft:SR  1,1.5;PR  -1,.5:LB2;PR  1,.5.SR  1.5,2;LB)" 

5644  END  IF 

5848  END  IF 
5852  IF  Opo-0  THEN 
5856  IF  Iun-0  THEN 

5860  PRINT  “SR  :.5,2;PU  PA  -12,38:01  0,l:LBh  /  (U/n;PR  -1,.5;SR  l,1.5iLB 

2:SR  1.5.2:PR  .5, .5: " 

5864  •  PRINT  “LB.:PR  .5,0;LBK)" 

5868  ELSE 

5872  PRINT  “SR  1.5,2;PU  PA  -12,28:01  0,l;LBh  /  ':8tu/hr;PR  - .  5 ,  .  5 ;  LB .  :  PR 

.5, .5: " 

5876  PRINT  "LBf t ; PR  -1 , .5; SR  1,1.5;LB2:SR  1.5,2;PR  .5,.5;LB.:PR  .5. .5; 
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LBF  i" 

5880  END  IF 

5884  END  IF 
5888  IF  ldl-0  THEN 

5892  PRINT  "SR  1.5,2:PU  PA  50,-16  CP" : -LEN< XlabelS )/2 ; “0; L8" ; XlabelS; 

5896  PRINT  "PA  -14,50  CP  0, " : -LEN( YlabelS  )/2*5/6 ; "DI  0, 1 ; LB" ; YlabelS; 

5900  PRINT  "CP  0,0  DI" 

5904  END  IF 
5908  Ipn-0 
5912  Repeat : ! 

5916  Xll-l.E+6 
5920  Xul  — l.E+6 
5924  Icn-0 
5928  BEEP 
5932  Ok-1 

5936  INPUT  "WANT  TO  PLOT  DATA  FROM  A  FILE  < 1-Y(DEFAULT ),0-N >?" ,0k 
5940  IF  Ok-1  THEN 
5944  BEEP 

5948  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" ,Df i le*( 0 > 

5952  ASSIGN  9File  TO  OfileS<0) 

5956  BEEP 
5960  Npaira>20 

5964  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED(DEFAULT-20 >" ,Npair9 
5968  BEEP 
5972  Itn»Itn+l 

5976  INPUT  "ENTER  TUBE  NUMBER  <1,  2,  3,  4.  OR  5)",Itn 

5980  BEEP- 

5984  PRINTER  IS  1 

5988  INPUT  "WANT  DEFAULT  SYMBOLS?  (YES-0  (DEFAULT  ),N0-1 )" ,Synb 

5992  Syii-Itn+2 

5996  IF  Synb-0  THEN 

6000  GOTO  6036 

6004  END  IF 

6008  PRINT  USING  '4X,  "Select  a  aynbol:"* 

6012  PRINT  USING  “6X,"1  Star  2  Plus  sion*” 

6016  PRINT  USING  “6X,""3  Circle  4  Square*" 

6020  PRINT  USING  "6X,"*5  Ronbus*" 

6024  PRINT  USING  ■GX,"6  Right-aide-up  triangle"" 

6028  PRINT  USING  "GX,""7  Up-aide-doun  triangle"" 

6032  INPUT  Syn 
6036  PRINTER  IS  705 
6040  PRINT  "PU  01* 

6044  IF  Syn-1  THEN  PRINT  "SM*" 

6048  IF  Syn-2  THEN  PRINT  "SM+* 

6052  IF  Syoi-3  THEN  PRINT  "SMo" 

G0SG  FOR  I-l  TO  Npaira 
6060  ENTER  0File;Xy(*  ) 

6064  IF  Opo-0  THEN 
6068  Ya-Xy(Itn-l) 
f072  Xa-Xy( 11  +  Itn  ) 

6077  END  IF 

6080  IF  Opo-1  THEN 

6084  Ya-Xy(5+Itn) 

6088  Xa-Xy( 11+Itn) 

6092  END  IF 

6096  IF‘Opo-2  THEN 

6100  Ya-Xy(5+Itn) 

6104  Xa-Xy(Itn-l> 

6108  END  IF 

6112  IF  Xa<Xll  THEN  Xll-Xa 
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6US  IF  Xa>Xul  THEN  Xul-Xa 

SI 20  IF  Iun-1  THEN 

6124  IF  Opo<2  THEN  Xa-Xa*1.8 

6128  IF  Opo>0  THEN  Ya-¥a«.17Gl 

6132  IF  Opo-0  THEN  Ya-Ya».317 

6136  IF  Opo-2  THEN  Xa-Xa».317 

6140  END  IF 

6144  X»L6T<Xa/Xnin)*Sfx 
6148  Y-L6T( Ya/Ynin>»Sfy 
6152  Kj-0 

6156  CALL  Synb(X,Y,Syn,Icl.KJ > 

6160  GOTO  6212 

6164  IF  Syn>3  THEN  PRINT  *SM“ 

6168  IF  Syn<4  THEN  PRINT  *SR  1. 4,2.4- 

6172  IF  Icl-0  THEN 

6176  PRINT  -PA-.X.Y,*- 

6180  ELSE 

6184  PRINT  -PA“,X,Y,*PD“ 

6188  END  IF 

6192  IF  SyM>3  THEN  PRINT  "SR  1. 2,1.6’ 

6196  IF  Syn-4  THEN  PRINT  -UC2,4,99,0,-8,-4,0,0,8.4,0i “ 
6200  IF  Syn-S  THEN  PRINT  •UC3, 0,99, -3, -6, -3, 6, 3, 6, 3, -6; " 
6204  IF  Syn-6  THEN  PRINT  ‘000,5.3,99,3,-8,-6,0,3.8! " 

6208  IF  SyKi-7  THEN  PRINT  *UC0.-5.3,99,-3.8,G.0,-3,-8: * 

6212  NEXT  I 

6216  PRINT  -PU- 

6220  BEEP 

6224  I lab-1 

6228  INPUT  -WANT  TO  LABEL?  ( 1-Y(DEFAULT >,8-N >• , Ilab 
6232  IF  I lab-1  THEN 
6236  PRINT  ■SP0iSP2' 

6240  BEEP 

6244  IF  Klab-e  THEN 

6248  X  lab-65 

6252  Y lab-85 

6256  INPUT  -ENTER  INITIAL  X,Y  LOCATIONS’ .Xlab.Ylab 

6260  Xtt-Xlab-5 

6264  Ytt-Ylab+8 

6268  PRINT  ’SR  l.l.S* 

6272  PRINT  ’SM-.PA"  ,Xtt  ,Ytt  .’LB  Tube  X  File" 

6276  Ytt-Ytt-3 

6280  PRINT  -PA’,Xtt,Ytt,"LB  No  Oil  Nane’ 

6284  IF  Syi-i-1  THEN  PRINT  ’SH-’ 

6288  IF  Syfi-2  THEN  PRINT  -SM+“ 

6292  IF  Syn-3  THEN  PRINT  -SHo" 

6296  Klab-1 

6300  END  IF 

6304  Kj-1 

6308  CALL  Synb(Xlab,Ylab,Syn, Icl.Kj ) 

6312  PRINT  -SR  l.I.SiSM’ 

6316  IF  Syn‘4  THEN  PRINT  ’PR  2,0” 

6320  PRINT  -PR  2,0.LB- ; Itn:  ” 

6324  BEEP 

6328  INPUT  -ENTER  BOP<0-OEFAULT >’ ,Bop 

6332  IF  Bop<10  THEN  PRINT  ‘PR  3,0:LB- :Bop; 

6336  ,  IF  Bop>9  THEN  PRINT  -PR  1.5,0iLB- iBop:  ” 

6340  PRINT  -PR  2 ,0iLB’ :0f ileS<0 >i * ’ 

6344  PRINT  -SP0:SPliSR  1.5,2’ 

6348  Ylab-Ylab-5 

6352  END  IF 
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6350 

6364 

6368 

6372 

6376 

6380 

6384 

6388 

6392 

5396 

6400 

6404 

6408 

6412 

6416 

6420 

6424 

6428 

6432 

6436 

6440 

6444 

6448 

6452 

6456 

6460 

6464 

6468 

6472 

6476 

6480 

6484 

6488 

6492 

6496 

6500 

6504 

6508 

6512 

6516 

6520 

6524 

6528 

6532 

6536 

6540 

6544 

6548 

6552 

6556 

6560 

6564 

6568 

6572 

6576 

6580 

6584 

6588 

6592 


ASSIGN  TO  • 

Xll-Xll/1,2 
Xul-Xul»1.2 
I  GOTO  8040 
END  IF 

PRINT  "PU  SM‘ 

BEEP 

Go_on“l 

INPUT  “WANT  TO  PLOT  A  POLYNOMIAL  <  1»Y(  DEFAULT  >  ,0-N  6o_on 
IF  Go_on«l  THEN 
BEEP 

PRINTER  IS  1 

INPUT  "WANT  DEFAULT  LINE  TYPE?  (YES-0  < DEFAULT ), NO-l) " .Ln 
Ipn-Itn 
IF  Ln-0  THEN 
GOTO  6448 
END  IF 

PRINT  USING  “4X,“  "Select  line  type:""" 

PRINT  USING  "6X,""0  Solid  line""" 

PRINT  USING  "GX,-"1  Dashed""" 

PRINT  USING  "GX.""2.,,5  Longer  line  -  dash""" 

INPUT  Ipn 
PRINTER  IS  705 
BEEP 
IIog-1 

INPUT  "SELECT  < 0-LIN,l-LOS< DEFAULT >)", Hog 
Iprnt-1 

CALL  Poly<Df ile$< *  ),Npairs,Itn) 

FOR  Xx-0  TO  Cx  STEP  Cx/200 
Xa-Xi>nn*10"Xx 

IF  Xa<Xll  OR  Xa>Xul  THEN  6572 

Icn-Icn+1 

Pu-0 

IF  Ipn-1  THEN  Idf-Icn  MOD  2 
IF  Ipn-2  THEN  Idf-Icn  MOD  4 
IF  Ipn-3  THEN  Idf-Icn  MOD  8 

IF  Ipn-4  THEN  Idf-Icn  MOD  16 

IF  Ipn-5  THEN  Idf-Icn  MOD  28 

IF  Idf-1  THEN  Pu-1 

Ya-FNPoly(Xa) 

IF  Ya^Ynin  THEN  6572 
IF  Iun-1  THEN 

IF  0po<2  THEN  Xa-Xa*1.8 
IF  Opo>0  THEN  Ya-Ya*.17Bl 
IF  Opo-0  THEN  Ya-Ya».317 
IF  Opo-2  THEN  Xa-Xa*.317 
END  IF 

Y-LGT( Va/ Ynin  )«Sfy 
X-LGT Xa/Xnin  )»Sf.-. 

IF  yv0  THEN  Y-0 

IF  Y:100  THEN  GOTO  6572 

IF  Pu-0  THEN  PRINT  "PA" ,X , Y, "PD" 

IF  Pu-1  THEN  PRINT  "PA" , X . Y . "PU" 

NEXT  X.-. 

.  PRINT  “PU" 

END  IF 
BEEP 

INPUT  "WANT  TO  QUIT  ( 1-Y .0-Ni DEFAULT  )>", Iqt 
IF  Iqt-1  THEN  6600 
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6596  SOTO  5916 

6600  PRINT  "PU  PA  0,0  SP0“ 

6604  SUBENO 
6608! 

6612  'XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
6616! 

6620  SUB  Synb(X,Y,Syn,Icl,Kj > 

6624  IF  Sypi>3  THEN  PRINT  'SM“ 

6628  IF  Syn<4  THEN  PRINT  “SR  1.4, 2. 4“ 

6632  Vad»0 

6636  IF  Kj-1  THEN  Yad-.8 

6640  IF  Icl-0  THEN 

6644  PRINT  “PA" , X , Y+Yad , " “ 

6648  ELSE 

6652  PRINT  “PA” ,X , Y+Yad, “PD“ 

6656  END  IF 

6660  IF  Syn>3  THEN  PRINT  “SR  1,2, 1.6“ 

6664  IF  Syn-4  THEN  PRINT  “UC2,4.99.0,-8,-4.0.0,8.4.0: “ 

6668  IF  Syn-S  THEN  PRINT  "UC3, 0.99, -3, -6, -3, 6. 3, 6, 3, -6; “ 

6672  IF  Syn-B  THEN  PRINT  “UC0, 5.3, 99, 3, -8. -6. 0.3. 8: “ 

6676  IF  Syin-7  THEN  PRINT  ’UC0, -5. 3, 99. -3, 8, 6. 0,-3, -8; “ 

6680  IF  Kj-1  THEN  PRINT  “SMiPR  0,-.8‘ 

6684  SUBENO 
6688! 

6692!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

66961 

6700  SUB  Fixup 
6704!  FILE:  FIXUP 
6708! 

6712  DIM  Enf(34  ),Anp( 11  ),Volt<4 >.Ldtc(4) 

6716  BEEP 

6720  INPUT  "OLD  FILE  TO  FIXUP"  ,02fiU« 

6724  ASSIGN  9File2  TO  02filet 
6728  OlfilaS-'TEST" 

6732  CREATE  BOAT  OlfiIeS,60 
6736  ASSIGN  eFilel  TO  OlfileS 

6740  ENTER  0File2i Nrun.DateS.LdtcC •  ), Itt .Bop.Nht .Natp.Nrt ,Corr 
6744  OUTPUT  0Fi lei • Nrun, OateC,Ldtc( •>, Itt .Bop.Nht ,Natp ,Nrt .Corn 
6748  FOR  I-l  TO  Nrun 

6752  ENTER  0Fi  le2  i  ToldS.Enf  (  •  ).  Volt(  •  )  ,At>ip(  *  ) 

6756  IF  I-l  THEN  6764 

6760  OUTPUT  9Fi lei i Bop , ToldS.Enf <•  ),Volt( •  ),Anp< • > 

6764  NEXT  I 

6768  ASSIGN  @File2  TO  • 

6772  ASSIGN  9Filel  TO  • 

6776!  RENAME  “TEST"  TO  02_file$ 

6730  BEEP  2000,. 2 
6784  BEEP  4000,.: 

6788  BEEP  4000,. 2 
6792  SUBEND 

6796' 

6800ixxxxxxxvxxxxxyxxxxxxxxxxxxxxxxxxxyxxxxxxKxxxxxxy)<xxxi'xxxxxxxxxxxxxxx 

6804' 

6308  SUB  Move 
68121  file  NAME:  MOUE 
6816' 

6820  DIM  A'’66  /  ,B(66  ),C(66  ),D(66  >,E(G6  ),F-:66  ),G(66  ),H(  66  ) .  J'  66  ),K(66  i.H  66  ),M(6B 
) 

6824  DIM  N(66,',Enf'34),Uolt<2),Anp(ll  .,Ldtc(4> 

6828  BEEP 
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6852  INPUT  "OLD  FILE  TO  MOUE" ,D2_f » le* 

6836  ASSIGN  SFileZ  TO  D2_filet 

6840  ENTER  «Fi le2 : Nrun, Toidt,Ldtc( •), Itt , Bop. Nht .Natp ,Nrt .Corn 

6844  FOR  I-l  TO  Nrun 

6848  ENTER  @File2iTold$ 

6852  ENTER  0File2; Ad  ) ,B( 1 ) ,C( I  ),D< 1 ),E( I ) .F< I ) ,6( I ) ,H(I ) . J< I > ,K( I > ,L( I > ,M( 
I),N<I) 

6856  ENTER  9File2 i Enf ( • > ,Uolt( • > ,Amp< • > 

6860  NEXT  I 

6864  ASSIGN  eFile2  TO  • 

6868  BEEP 

6872  INPUT  “SHIFT  DISK  AND  HIT  CONTINUE", Ok 
6876  BEEP 

6880  INPUT  "INPUT  BOAT  SIZE". Size 
6884  CREATE  BOAT  D2_f ileS.Size 

6888  ASSIGN  @Filel  TO  D2_fiie$ 

6892  OUTPUT  9Fi lei i Nrun, Date*, Ldtc< •), Itt , Bop, Nht ,Natp ,Nrt ,Corr 

6896  FOR  I“1  TO  Nrun 

6900  OUTPUT  9Fi  lei -.Told* 

6904  OUTPUT  0Filel  i  A(  I  >,B(  I  ),C<  I  ),0(  I  >,E(  I  >,Fd  >,6(  I  ),H<  I  >,  J(  I  ),K(  I  >,Ld  ),M 
(I).Nd) 

6908  OUTPUT  9Fi lei • Enf ( * ) ,Uolt < •  ) , Anp( * > 

6912  NEXT  I 

6916  ASSIGN  0Filel  TO  • 

6920!  RENAME  “TEST"  TO  D2_file» 

6924  BEEP  2000,. 2 
6928  BEEP  4000,. 2 
6932  SEEP  4000,. 2 
6936  PRINT  “DATA  FILE  MOVED" 

6940  SUBENO 
6944! 

G948IXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6952* 

6956  Sub  Purg 
6960  BEEP 

6964  INPUT  "ENTER  FILE  NAME  TO  BE  DELETED" ,Fi le* 

6968  PURGE  File* 

6972  GOTO  6960 
6976  SUBENO 

6980! 

6984!XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

6988! 

6992  SUB  Conb 

6996!  FILE  NAME:  COMB 

7000! 

7004  DIM  Enf<34),Volt(2),Anp(ll),Ldtc<4> 

7008  BEEP 

7012  INPUT  "OLD  FILE  TO  FIXUP" .D2_file* 

7016  ASSIGN  9Fiie2  TO  D2_file* 

7020  Dl_f ile*-”TEST” 

7024  CREATE  BOAT  Dl_file*,30 
7028  ASSIGN  SFilel  TO  Dl_file$ 

7032  ENTER  SFi le2 ; Nrun,Oate*,Ldtc' •). Itt .Bop.Nht .Natp ,Nrt , Corn 

7036  IF  K“0  THEN  OUTPUT  @Fi lei i Nrun ,Oate* ,Ldtc( •  ) , 1 1 1 ,Bop ,Nht , Natp , Nrt  ~orr- 

7040  FOR  I-l  TO  Nrun 

7044  ENTER  @Fi  le2 :  Bop  ,  Told*  ,Enf(  •),  Volt  (•)  ,Apip(  •  ) 

7048  OUTPUT  0Fi lei : Bop , Told* , Enf( •) ,Uolt( •) ,Anp( •  ) 

7052  NEXT  I 

7056  ASSIGN  9File2  TO  • 

7060*  RENAME  "TEST"  TO  D2_file* 
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7064  BEEP  4000,. 2 
7068  BEEP 
7072  Oka-1 

7076  INPUT  “WANT  TO  rtOD  ANOTHER  FILE  < 1-Y,0-N( default >>?“, Oka 
7080  IF  Oka-1  THEN 
7084  K-1 

7088  BEEP 

7092  INPUT  -GIUE  NEW  FILE  NAME” ,Nf i le$ 

7096  ASSIGN  eFile2  TO  NfileS 

7100  GOTO  7032 

7104  END  IF 

7108  ASSIGN  9File2  TO  • 

7112  SU6EN0 
7116! 

71201 xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

7124! 

7128  SUB  Readplot 

7132  DIM  0dp(5),Htube<5),Theiab<5) 

7136  PRINTER  IS  701 

7140  INPUT  "ENTER  FILE  NAME'.Fxle* 

7144  INPUT  “ENTER  THE  NUMBER  OF  DATA  PAIRS’, Nrun 
7148  ASSIGN  eFilal  TO  File$ 

7152  FOR  I-l  TO  Nrun 

7156  ENTER  CFxlel i Odp< *  ) ,Htube( • ),Thetab( * > 

7160  PRINT  Qdp(*> 

7164  PRINT 

7168  PRINT  Htube<*) 

7172  PRINT 

7176  PRINT  Thetab<») 

7180  PRINT 

7184  PRINT 

7188  NEXT  I 
7192  SUBENO 


175 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 

1 .  Defense  Technical  Infonnation  Center  2 

Cameron  Stadtrn 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  California  ^3943-5000 

3 .  Professor  Paul  J.  Marto,  Code  MEMX  2 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

4.  Professor  Steven  Memory,  Code  MEMR  1 

Department  of  Mechanictd  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

5.  Department  Chairman,  Code  ME  1 

Departmoit  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  California  93943-5000 

6.  Naval  Engineering  Curricular  Officer,  Code  34  1 

Department  of  Mechanical  Engineering 

Naval  Postgraduate  School 
Monterey,  Cilifomia  93943-5004 

7.  Mr.  R.  Helmick,  Code  2722  1 

£)avid  Taylor  Research  Center 

Annapolis,  Maryland  21402-5067 

8.  Mr.  Bruce  G.  Unkel  1 

NAVSEA(Code56Y15) 

Department  of  the  Navy 
Washington  D.C.  20362-5101 

9.  LT.  Scott  V.  Chilman  1 

720  Hill  Rd. 

Philadelphia,  PA.  19128 


176 


